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Chapter 1 
Introduction 



'Modern physics', as it stands one century after Einstein's annus mirabilis, has had 
the privilege of giving scientific footing and empiric evidence to some insights that 
philosophers had long ago. Immanuel Kant already stated in his "Critique of Pure 
Reason" [1] that space and time are inherent concepts of our mind, not properties of 
the objects we perceive. These ideas always seemed too abstract to be taken seriously 
by non-professional philosophers. But nowadays students in special relativity and 
quantum mechanics, optical engineers, molecular biologists all at some point in their 
studies face the severe difficulties that a naive realistic world view (possibly even 
with a notion of universal time) has in the light of the past 100 years of science. We 
become painfully aware how much our thinking depends on concepts such as space, 
time and causality- prerequisites for our mind to work. They are given to us "a 
priori", and we can do nothing but to accept that they don't seem to be adequate 
in the context of modern physics. 

Nowadays, we're getting used to lasers, the global positioning system and ele- 
mentary particles such as electrons. But while we have gained a hundred years of 
experience in dealing with quantum mechanical and relativistic laws and "objects" , 
the basic insights of modern physics remain astounding. 

The three-loop quark form factor that this thesis deals with is a relativistic 
quantum mechanical amplitude. It describes a quantum mechanical probability 
amplitude (it is a complex number) and it is a real physical object in the same 
sense as electrons are (quantum mechanical) waves (not classical "spheres"). It is 
tempting to think of the quark form factor as a part of the classical probability of 
an electron scattering off a proton, rather than of a quantum mechanical amplitude. 
Classical intuition is comfortable, but it must be kept in mind that it is just that: 
comfort for a mind that is not adapted to the laws of nature it explores. 

Is it real? To the critical reader and also to non-experts it must be admitted that 
modern experimental particle physics has become extremely complex, and therefore 
measurements are more indirect. To measure the quark form factor, modern par- 
ticle physics experiments such as LEP, HERA, and very soon the eagerly awaited 
Large Hadron Collider (LHC) study various processes such as jet production, deep- 
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Figure 1.1: The tree-level and one- loop diagrams contributing to the quark form 
factor. 

inelastic scattering, Drell-Yan photo-production etc. A measurement of a physical 
quantity is a very complex task involving (just to give a rough and incomplete list) 
fund raising and constructing the facility, testing and calibrating the experiment, 
performing Monte-Carlo simulations of the collisions in the detector and finally 
comparing experimentally measured vs. theoretically expected outcomes of the ex- 
periment. Nevertheless, to particle physicists, the quark form factor must be real 
simply because it is a gauge invariant quantity in quantum chromo dynamics (QCD), 
and because the QCD-improved parton model has been tested so successfully in the 
perturbative regime. But such an experts answer is clearly neither exhaustive nor 
the only possible one. And since almost the entire rest of this thesis will strive to 
be professional and answer all questions about the form factor by putting a single 
number to it, I'd like to take this moment and pass the question on to someone 
with an entirely different viewpoint. Maybe, if physicists claim the form factor to 
be real, it should be real also for non-physicists? Maybe it isn't just a number? "If 
the quark form factor is real, how does it look like?"- Please turn the page and let 
an artist give her answer. 

Since it is a QCD quantity, the immediate context of the quark form factor 
is not so much new physics but rather a precise understanding of "known" physics 
needed to discern small new physics deviations from known physics uncertainties. In 
perturbation theory the quark form factor is expanded in small coupling constants, 
the largest of which being ag. The contributions to various orders in this expansion 
(called "tree-level", "one-loop", "two-loop") can be represented as Feynman graphs. 
To physicists, Feynman graphs are useful because the converse is also possible: Each 
Feynman graph corresponds to a quantifiable contribution to the form factor. The 
tree-level and one-loop [2] contributions shown in fig. 1.1 are text-book calculations, 
and the two-loop result in fig. 1.2 has been obtained more than 20 years ago [3]. The 
three-loop contribution is expected to be a correction of a few percent. Observing 
it at the LHC might thus be difficult, but possible. For the International Linear 
Collider (ILC), it clearly is an interesting object, helping to bring the theoretical 
prediction of several observables to percent-level precision. 
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Figure 1.2: Diagrams contributing to the two-loop quark form factor. The eUipsis 
denotes 13 further diagrams. 

In this thesis, only part of the calculation is attempted. Chapter 2 will use 

standard tools such as Qgraf [4] and FORM [5] to generate the amplitude as a 
linear combination of three-loop integrals. In chapter 3 we will then have to anal- 
yse these integrals closely. When first obtained, they are in a quite unstructured 
form and we have to devise a (partially) automated way to bring these integrals 
into a unique standard form. This is necessary because in chapter 4 we will solve a 
system of linear equations in those integrals in order to (dramatically) reduce the 
number of integrals that actually have to be evaluated. These systems of equa- 
tions would be unnecessarily large without standardisation of the integrals. In 
view of the fact that so far we have not managed to obtain the final result for 
the reduction of the three-loop form factor, chapter 4 will show why this task is 
so hard, despite the fact that a CPU cluster was at our disposal (http://www- 
theorie.physik.unizh.ch/ stadel/zbox/start7^zbox2). We will then introduce two im- 
mediate improvements (symbolic reduction of reducible topologies and loop-by-loop 
integration), and also report on improvements on the reduction algorithm we use. 




Is the three-loop QCD form factor real? If so, how does it look like? An artists impression, courtesy E. Schiirmann. 



Chapter 2 

The Quark Form Factor 



In this chapter we will introduce the three-loop QCD quark form factor. This is the 
object we ultimately wish to compute. We will also motivate why its computation 
should be worth the effort and where it could improve the theory prediction of 
physical observables. 

2.1 Introducing the Quark Form Factor 

In perturbative QCD, the quantum mechanical amplitude for the process 7* qq, 
M. (7* — > qq) is expanded in the strong coupling constant a^. It is then represented 
as a sum of Feynman graphs: 



The shaded blob represents the sum over all possible higher order corrections that 
are allowed by the Feynman rules of QCD. The matrix element M.^ (7* qq) has 
an incoming off-shell photon and an outgoing quark-antiquark pair. 

In a theory like QCD that we take to be CP-even, the only vectors that can 
carry the Lorentz index of the external photon are the Dirac matrices 7'* and the 
(outgoing) momenta and of the external quark and antiquark, respectively. 
So we may write [6] 



MM,-' ^ qq) ^ (_iee,) u^{p,) [A^ + Bq^ + C K - p^)] v^' (p^) , (2.1) 



where q = Pi + P2 and Cg is -1-2/3 for up- and —1/3 for down-type quarks, u and v 
are the spinors for the quark and the antiquark, respectively, and they depend on 




q 



7 



* 



q 



9 



10 CHAPTER 2. THE QUARK FORM FACTOR 

the particle's spin and momentum. If we considered the process Z* — > qq, couphngs 
proportional to 7'^75, q^^5 and ia^^qv would also be present. 

Using the Ward-Identity q^Ai'^ = and the Dirac equations of motion, we 
find that B = 0. Finally, in massless theories with vector-like gauge bosons (like 
QED and QCD), helicity is conserved, which leads to C = 0. Taken together, the 
massless vertex-corrections Ai^ (7* —>■ qq) is given by a single form factor T = A. 
Using Lorentz covariance of Al^ and 7'* in (2.1), it must be a Lorentz scalar and 
thus it may only depend on kinematical invariants. Since the only non-vanishing 
kinematical invariant is g^, we have 

A^r' (7* - m) = (-^ee,) T {q') u^ip,) 7X'(p2) • (2.2) 
The form factor may now be extracted from (2.2) via 

UfUc 

Q 1 s,s 

Eq. (2.3) computes the quantity we want, {q^), out of M.'^'^^' (7* — > qq), which is a 
sum of three-loop Feynman graphs. In sec. 3.1 we will describe their generation and 
how to evaluate the sum over spins in (2.3). The rest of chapter 3 will be concerned 
with finding a compact standard representation for the three-loop integrals that we 
obtain after this summation over spins. This is preparatory work necessary for the 
reduction to master integrals that will be described in chapter 4. 

2.2 Physics Motivation for the Quark Form Fac- 
tor 

Having defined the quark form factor T {q^), we next have a look at its physical rel- 
evance. It constitutes part of the calculation for various different physical processes, 
on some of which we shall comment below. Apart from this its infrared behaviour is 
of special interest since it may be used to determine various (process independent) 
quantities in the context of resummation. 

The three-loop contribution, being formally O {aa^) , is expected to contribute 
a correction of a few percent and thus in many cases it will probably increase the 
precision of the theory prediction to the percent level and possibly beyond. The 
corresponding accuracy will be reached in experiments both at the LHC and the 
ILC. 

2.2.1 The Form Factor in e+e" Collisions 

The three-loop QCD form factor directly enters the matrix element for e+e~ — > qq. 
It may be written down using eq. (2.2), where the form factor acts as an effective 
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7-g-g-vertex and thus includes the higher order QCD corrections. This effective 
vertex is denoted by a blob in fig. 2.1. Clearly, if three- loop QCD corrections are 
included in this way, QED corrections cannot be neglected, and they may not be of 
the form of fig. 2.1. For example, an additional photon may be exchanged between 
the external electron and quark line. 

Now the cross section for e'^e~ qq is not infrared finite at higher orders in the 
perturbative expansion, and neither is the form factor. Thus so far the QCD form 
factor will only summarise an important contribution to physical processes. Other 
contributions come from radiative corrections with additional external gluons and 
photons. Physical observables are e.g. the cross sections for e+e~ — > n jets and all 
their differential distributions (event shapes, angular and pr-distributions etc.). At 
present they are known up to (at most) NNLO, and there the contribution from the 
two-loop vertex correction has certainly been one of the easier parts of the calcu- 
lation. The N^LO contribution, to which the three-loop form factor would belong, 
is clearly not yet within reach of a theoretical computation. Wc may still view the 
three-loop contribution as a first step in this direction. Knowing the theory value 
of the form factor to three loops should eventually lead to an improved determi- 
nation of the strong coupling constant and its running at the ILC. New physics 
can then be envisaged once the standard model parameters are known precisely 
enough. For example, the present error on represents the dominant uncertainty 
on the prediction of the scale for grand unification of the strong and electro-weak 
interactions [7]. 




Figure 2.1: The form factor enters e+e — > 2jets. 



2.2.2 The Form Factor in Deep-Inelastic scattering 

Deep-inelastic scattering (e.g. at HERA, Hamburg) prominently involves the ex- 
change of a hard photon between the incident lepton and the target hadron. At 
parton level, the photon is incident to a quark line and so the quark form factor 
enters directly here (see fig. 2.2 for an illustration). Deep-inelastic scattering experi- 
ments usually measure structure functions which then allow to determine parton dis- 
tribution functions. Together with Drell-Yan lepton pair production deep-inelastic 
scattering processes are among the best suited ones for probing the short-distance 
structure of hadrons. The uncertainty from higher order perturbative corrections 
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often dominates and so higher order calculations like this three-loop correction are 
needed. 




q 



Figure 2.2: In deep-inelastic scattering, the exchanged photon probes the target 
nucleon structure via the form factor. 



2.2.3 The Form Factor in Drell-Yan Lepton Pair Produc- 
tion. 

Data on Drell-Yan lepton pair production, pp i^fr (or pp — > £^£^) will reach 
unprecedented precision at the LHC. It will clearly be possible to see corrections 
of the order of one percent or even in the per-mille range for certain observables 
such as the W- or Z-boson mass and sin^vi'- Fig- 2.3 shows how the q-q-^ form 
factor enters the Drell-Yan cross section. To make quantitative predictions also for 
differential distributions in the Drell-Yan process, the parton distributions for both 
incident ( ant i-) quarks are needed: 

^ - E JJ^^ disfalA iU, /^) ^/VB i^B, /.) (2.4) 

Here, fa/A denotes the parton distribution function and can be interpreted as the 
probability to find a parton of type a in the hadron of type A carrying a fraction C,a 
to + d^A of the hadron's momentum. They are universal in the sense that they 
are independent of the hard scattering process. They are best extracted from deep- 
inelastic scattering experiments. Due to the uncertainties in the parton distribution 
functions, Os-measurements will not reach the precision of LEP or ILC experiments. 
Rather, as stated before, the form factor will be an ingredient of "known physics" 
that will be needed when searching for "new physics" . 



2.2.4 The Form Factor in the Context of Resummation 

Form factors also play an important role in a different way: Using factorisation and 
gauge invariance it can be shown that infrared singularities of the form In" {q^ / 1^^) 
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Figure 2.3: The Form Factor enters the Drell-Yan cross-section when a virtual pho- 
ton is exchanged. 

can be resummed. They exponentiate according to [8] 




By matching the general exponentiated expression for the form factor onto a fixed 
order calculation, rcsummation coefficients in the exponent may be determined. 
These in turn allow -after reexpansion of the exponential- to determine subleading 
logarithms in the coefficient functions for processes such as deep-inelastic scattering 
or Drell-Yan lepton pair production. 

2.3 Current Status of QCD and Electroweak Cor- 
rections to Form Factors 

2.3.1 QED Form Factors 

Vertex corrections are prominent already in typical first year quantum field theory 
courses when the O (a) correction to the anomalous magnetic moment of the elec- 
tron [9] is studied. The higher order contributions to the QED vertex can be used 
to determine the fine structure constant a. These calculations have been pushed to 
O (a^) and turned this observable famous [10, 11] for yielding one of the most pre- 
cise predictions in physics. The muon anomalous magnetic moment is particularly 
interesting for finding new physics effects [12]. Calculations for massive two- loop 
vertex corrections have also been performed already [13]. 
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2.3.2 QCD Form Factors 

The one- loop QCD vertex correction was derived in 1949 by Schwinger [2]. The 
massless two-loop contribution became available by the work of Gonsalves [3] in 
1983, and also by Kramer and Lampe [14] and van Neerven et al.[15]. Massive 
external quarks were also considered at two loops. As was already mentioned in 
sec. 2.1, there are four (CP even) form factors in this case, two vector and two axial 
vector form factors. They were calculated in [16, 17, 18]. 

Quite recently, in two by now famous papers Moch, Vermaseren and Vogt com- 
puted the three-loop QCD sphtting functions [19, 20], from which they derived the 
divergent part (the 1/e poles in dimensional regularisation) of the three-loop quark 
form factor [21]. 

2.3.3 Elect roweak Form Factors 

Due to the additional scales of the W- and Z-boson masses, electroweak correc- 
tions to the form factor are much more difficult. One-loop results for observables in 
e'^e~ ff (where / is either a lepton or a quark) have already been available for 
LEP [22]. Full two- loop results do not seem to be in reach yet, though progress is 
made [23]. It was then reahzed that leading and sub- leading "Sudakov"-type loga- 
rithms, \nq'^/m^/z contribute corrections of up to 10% [24, 25] at — 1 TeV. Thus 
aiming at predictions for the ILC the community has focused on resumming (leading 
and next-to-leading) logarithms of over Mw,z for various processes (see [26] and 
references therein). 



Chapter 3 

Generation and Standardisation 



We now concentrate on the three-loop QCD contribution to the quark form factor. 
As a first step, section 3.1 will explain how to obtain the form factor as a hnear 
combination of three-loop integrals. Once obtained, these Feynman integrals will 
still be quite unstructured. Sections 3.2, 3.3 and 3.4 then describe how the following 
goals are met: 

• To find a compact representation for all three-loop vertex integrals, 

• To identify those topologies that vanish trivially, and 

• To exploit all (obvious) symmetries to make this representation as unique as 
possible. 

In this way, we will arrive at a unique representation of any three-loop integral that 
we will have to compute. The actual computation is attempted in chapter 4. 

3.1 Generating the Three-Loop Quark Form Fac- 
tor 

We would like to obtain the form factor using eq. (2.3) by providing the matrix 
element M.'^'^'^' (7* — > qq) exphcitly. 

It is well known that Feynman diagrams exhibit factorial growth in external legs 
and also in loops. We will soon see that for the three-loop form factor, one faces 
-depending on how diagrams are counted- several hundred diagrams. Thus it is 
clearly necessary to automatise the generation and evaluation of these diagrams. A 
number of tools have been written to this end, e.g. FeynArts [27] and Qgraf [4]; in 
appendix A we give further details on how Qgraf and FORM [5] were used for the 
generation of the quark form factor. 

Qgraf essentially requires the input of a model file, a style file, and the file 
qgraf . dat. The model file describes which types of vertices and propagators appear 
in the diagrams to be generated. Qgraf output is user-defined according to the 
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style file. When Qgraf is run, it reads the process definition given in qgraf . dat 
and then generates information on each possible diagram. The style file then defines 
how this information is written to the output, a plain text file. The style file may 
be tailored such that the output can be read into FORM. 

When Qgraf is run it generates the three-loop diagrams for 7* —>■ qq. The output 
of Qgraf contained 244 diagrams. 

Prom here, the Feynman amplitude may be constructed in FORM. The projec- 
tion onto the form factor, eq. (2.3), involves summing over spins and colour factors. 
They lead to traces of Dirac and colour matrices, respectively. FORM has built-in 
capabilities for evaluating Dirac traces, and colour traces are easily evaluated using 

fabc ^ zlxr , (3.1) 

Tr 

with Tr — 1/2. The Dirac traces produce scalar products involving loop and external 
momenta. The propagators involve the loop momenta in a way dictated by Qgraf. 
For example, the contribution of the diagram discussed in more detail in appendix A, 
fig. 3.2, will involve integrals of the form 



h ■ k2, ih ■ qf ,etc. 

kfklkl {{ki +pifY {ki + qf {k, -k2- ks+p,f ' 

In order for this output to be useful, we will next exploit the freedom to shift the 
loop momenta, in such a way as to bring the propagators into a unique form. The 
next chapter will be devoted to this task. 



3.2 Representing Integrals in a Compact Way 

In this section our goal is to set up a compact notation that is capable of representing 
any possible three-loop integral that we might encounter in our reduction. As a 
first step, we will have a look at the possible scalar products that appear in the 
numerators and denominators of our integrals. 

The three- loop vertex- integrals have denominators consisting of propagators, and 
numerators that are scalar products stemming e.g. from the Dirac traces of the 
quark line. Numerators and denominators consist of scalar products involving the 
loop momenta ki, k2, ks and external momenta pi and pi2 — Pi +P2- Out of these, 
13 different nonvanishing scalar products can be built: 

ki ■ ki k2- k2 ks • ks 

ki ■k2 ki ■ ks k2 ■ ks 

ki -pi /C2 ■ Pi ks ■ pi (3.4) 

ki ■ Pl2 k2 ■ Pl2 ks ■ Pi2 

Pl2 ■ Pl2 ■ 
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These scalar products span a (real) vector space. The set of possible inverse 
propagators that might occur in our three-loop integrals are elements in it. So 
alternatively, we may choose as our basis any set of inverse propagators, e.g.: 

1.2 1.2 1.2 

1^2 "'S 

{kl - k2f {kl - k^f {k2 - k^f 

{ki-pif {ki-pi2f (3.5) 

{k2-pif {k2-pi2f 

{ks-pif {ks-puf 



together with p^2- The basis transformation from (3.4) to (3.5) reads 



kl 


■kl 


^1 




k2 


■k2 






ks 


■ks 






kl 


■k2 


-[kl- 


{ki-k2f + kl]/2 


kl 




= [kl- 


{ki-ksf + kl]/2 


k2 


■ks 


= [kl - 


{k2-ksf + kl]/2 


kl 


■Pi 


= [kl- 


{ki-pif]/2 


kl 


■P2 


= [{(kl 


-Pl)') -(fci-Pi2)'+pL]/2 


k2 


■Pl 


= [kl - 


{k2-pif]/2 


k2 


■P2 


= [ik2 - 


-Plf - (A;2 -P12)' +Pl2]/2 


ks 


■Pl 


- [kl - 


{ks-pif]l2 


h 


■P2 


= [ik3 - 


- Plf - {ks - Puf + P12] / 2 ■ 


Pi 


■P2 


=Pl2/2 





(3.6) 



Now a Feynman integral is determined if the propagators and scalar products 
are given as a linear combination of (3.4). Alternatively, we may use (3.5) together 
with p^2- In general, each propagator of an arbitrary graph will be given as a linear 
combination of the propagators in (3.5) and pf2- 

Still aiming towards a compact representation of the three-loop vertex integrals, 
we will now argue the following: For integrals occurring in planar diagrams, it is 
possible to shift loop momenta such that all occurring propagators are members of 
the basis set (3.5), rather than just a linear combination of them. For example, 
the scalar integral (the one without any scalar products in the numerator) for the 
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q 



Figure 3.1: A sample planar diagram, called i. It is the only t = 9 planar master 
integral that appears in the reduction of the three-loop quark form factor. 



diagram in fig. 3.1 can be represented as 



So the integral is now given by a compact expression indicating the 12 exponents 
of the (inverse) propagators (3.5). If the 'propagator' l/pi2 should ever appear, we 
will simply factor it out of the integral and write it down explicitly. Its exponent 
could even be inferred by dimensional analysis, since it constitutes the only scale 
left after the loop integration is performed. 

The notation introduced in (3.7) includes the possibility of having the same 
propagator appearing several times, as is the case in the diagram of fig. 3.2, and 
in many integrals generated in integration- by-parts identities (see below, chap. 4). 
Furthermore, if the numerator of the integral contains scalar products, they may be 
rewritten as inverse propagators using (3.6). 

Thus this notation is capable of expressing any three-loop integral (including 
scalar products) if its propagators are already shifted such that they are members 




1 



X 



{h - pif {k2 - qf {h - Pif ih - qf 

= r(-l, -1, 0, -1, -1, -1, -1, 0, 0, -1, -1, -1) . (3.7) 
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Figure 3.2: A diagram whose topology inherently consists of two propagators with 
the same momentum flow. 

of a standard set. As a more general example for the same diagram fig. 3.1, we have 

/ d%d%d% , 

J (klf kl (k^ - k^f {{k^ - ksf) (k, - k^f 

1 

(^1 - Pi) (^2 - q) {ks - Pi) {ks - q) 
= ^ ( - ri(-3, -1, 0, -2, -1, -1, -1, 0, 0, 0, -1, -1) 

+ ri(-3, -1, 0, -2, -1, -1, -1, 0, 1, -1, -1, -1) (3.8) 
+ pL^i(-3, -1, 0, -2, -1, -1, -1, 0, 0, -1, -1, -1)) , 

where the subscript 1 in 7i shall specify that the propagators are taken with respect 
to the propagators in (3.5). We will refer to a (standard) set of propagators such 
as (3.5) as an auxiliary topology. A topology is specified by indicating a subset of 
distinct propagators of the auxiliary topology. 

3.2.1 Auxiliary Topologies 

One important virtue of the notation introduced in the previous section lies in the 
fact that usually many different Feynman graphs can be represented in one auxiliary 
topology. For example it turns out that all planar^ diagrams share the auxiliary 

^The notions "planar" or "nonplanar" are properties defined in graph theory. We use them 
without much formahty. Planar diagrams arc diagrams that can be drawn on the plane in such 
a way that propagators only intersect at vertices. Nonplanar graphs are graphs where this is not 
possible. 



20 



CHAPTER 3. GENERATION AND STANDARDISATION 




Figure 3.3: A sample non-planar diagram belonging to the auxiliary topology num- 
ber 2. 



topology already mentioned earlier: 

1.2 1.2 1.2 

tvi Atg 

{ki - {ki - ksf {k2 - ksf 

{ki-pif {ki-puf (3.9) 

ik2-pif {k2-pi2f 

{ks-pif {ks-puf- 

In general however, a single auxiliary topology will not be enough. In particular, 
non- planar diagrams do not fit into the propagator set (3.9). As an example, the 
diagram in fig. 3.3 has the propagators 



kl, kl, (ki - kif , {ki - k^f , {ki -k2- k^f 
(ki - pif , (kl - qf , {k2 - pif , {h - qf ■ 

Clearly, the inverse propagator [ki — k2 — k^)^ in (3.10) is not one of (3.9). So 
it is necessary to define a new auxiliary topology that contains this propagator and 
again has 12 hnearly independent propagators. We choose 

h.2 h.2 h.2 

tv^ tV2 'vg 

{kl - k2f {kl - ksf {kl - k2- ksf 

{ki-pif {ki-puf (3.11) 

{k2-Pif {k2-Pi2f 

{h-Pif {h-Puf- 

Finally, a third auxiliary topology is needed to represent the propagators of diagrams 
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Q 



q 

Figure 3.4: A sample non-planar diagram belonging to the auxiliary topology num- 
ber 3. 



ID 


Topology 1 


Topology 2 


Topology 3 


1 


h 


ki 


ki 


2 




k2 


k2 


4 






ks 


8 


ki - k2 


ki - k2 


ki - k2 


16 


ki - ks 


ki - ks 


ki - kz 


32 


k2 - ks 


k2 — k^ 


ki-k2- kz 


64 


ki -pi 


ki- ks- p2 


ki -pi 


128 


ki - Pi2 


ki - Pi2 


ki - Pi2 


256 


k2-Pi 


k2-Pi 


k2 -Pi 


512 


k2 - Pl2 


ki-k2- P2 


k2 - Pl2 


1024 


k^-pi 


ks-pi 


ks-pi 


2048 


k3 - Pl2 


ks - Pl2 


ks - Pl2 




Table 3.1: The auxihary topologies needed to represent all integrals appearing in the 
three-loop form factor. Topology 1 corresponds to planar diagrams. The numbers 
on the left will be convenient to determine the topology ID (see sec. 3.4 below). 
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like the one in fig. 3.4. Two new propagators are needed: 1/ {ki — ks + Pi f and 
l/{ki-k2 + P2) . We may choose the auxihary topology 



h.2 h.2 h.2 

tvi tv2 ^3 

(fci - fcs)^ ih - kzf {k2 - ksf 

(h-ks-p2f (ki-puf ^ (3.12) 

{k2-pif {ki-k2-p2f 

{h-Pif {h-Puf- 



We find that all integrals appearing in the three-loop quark form factor can be 
represented in one of the auxiliary topologies (3.9), (3.11) and (3.12). Tab. 3.1 sum- 
marises the auxiliary topologies that are needed to represent all integrals appearing 
in the thrcc-loop quark form factor. 

In what follows, we will concentrate on the case of planar diagrams. An omitted 
index 1 on integrals T = Ti will be understood. For several reasons, the reduction 
of planar topologies is the most important and most difficult one: Firstly, 202 out 
of the 244 diagrams will turn out to be planar below. Furthermore, the reduction 
algorithm discussed in chapter 4 will show that integrals of a given topology can be 
reduced to integrals of so-called subtopologies. They have one type of propagator less 
than the parent topology. Since for the quark form factor there are no non-planar 
topologies with less than seven propagators, the non-planar topologies will reduce 
to planar ones at the level of six-propagator topologies. In addition quite a few 8- 
or 7-propagator subtopologies of non-planar topologies are planar. 

3.2.2 Casting Propagators into a Standard Form 

The propagators of the Qgraf output differ from those in tab. 3.1 by shifts in the 
loop momenta. For example, the Qgraf output for the diagram in fig. 3.2 contains 
integrals such as the ones given in (3.3). By inspection, one finds that after applying 
the shift 



all propagator momenta are members of the auxiliary topology 1 in tab. 3.1, up to 
a sign. Thus in order to bring the integrals of the Qgraf output into the compact 
form defined in the last section, the task is to find this shift and apply it to the 
Qgraf output; to rewrite the scalar products as inverse propagators and cast all 
integrals into the form of one of the auxiliary topologies in tab. 3.1. 

How to find these shifts in the loop momenta? Formally we are given t propaga- 
tors that arc expressed as linear combinations of k = {ki, k2, k-i,pi,p2}; i.e. the set 
of propagator momenta can be written as A.k where A is a t x 5-matrix. Similarly 
the target set of propagator momenta can be written as B.k'. We are looking for a 



ki 
k2 

h 
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change of basis k' — S.k such that A.k — B.k' up to a permutation P, and up to a 
sign in each component. I.e. given A and B, find a permutation P and a change 
of basis S such that A.k = P.B.S.k up to a sign in each component. A brute force 
algorithm to find S and P can be set up hke this: 



1. Choose any P. 

2. Set up the system A.k = P.B.S.k. 

3. Pick the first equation out of the system that contains ki and extract the loop 
momentum shift, ki — Suk^. 



4. Apply this shift to the system, A.k 



and check whether 



= P.B.S.k 

it might still have a solution. I.e. there should be no equation that does not 
depend on unprimed loop momenta and that is not a member of B.k'. 



5. If a solution is still possible, repeat step 3 and 4 with loop momentum k2 
(instead of /ci), then k^. If not, choose the next P and go to step 2. 

6. Upon arriving here, ki — > Syk^ is a solution. 

The problem is the permutation P- because naively trying all 12! ~ 5 • 10^ possi- 
bilities is out of question. But the check in step 4 prevents a lot of permutations 
of being fully tested. The average computation time for a given diagram and with 
B according to the set (3.9) was ~ 15 seconds on a 1.7GHz Pentium processor. 
Diagrams where no solution could be found were all non-planar. 

For planar graphs a different (and faster) approach is worth pointing out: There, 
a unique shortest path between the external legs exists. Then, the following algo- 
rithm can be used: 

1. Find the shortest paths between any two of the three external vertices. 

2. Out of these three paths, choose the longest one (or any one of the longest 
ones) . 

3. Prescribe the following loop momenta on the first two legs: 

• For the j-q leg: ki, k2] 

• For the 7-g leg: ki - pu, k2 - Pu] 

• For the q-q leg: ki — pi, k2 — pi. 



4. 



Shift the third loop momentum to fit a chosen propagator into any of the 5 
/ca-dependent propagators. Try all possibilities if necessary. 
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In the above procedure, steps 1 to 3 fix two of the three loop momenta. This is 
always possible because there is no nonvanishing three-loop vertex graph whose 
longest shortest path between external vertices contains just one propagator^. For 
the third loop momentum, at most five shifts have to be considered (since in the 
planar auxiliary topology each loop momentum occurs in five propagators). The 
only danger here might be snail (or tadpole'^) diagrams; but Qgraf has an option to 
switch off the generation of those diagrams. 

3.3 Eliminating Vanishing Topologies 

In this section we will identify those loop integrals that evaluate to zero. The re- 
duction using integration-by-parts identities will prove difficult and computationally 
expensive, so topologies that are zero should not be treated in this way, nor should 
they be kept in reduction equations of other topologies. 
In dimensional regularisation, integrals of the form 

vanish. They do not depend on any scale and therefore must be zero, simply because 
applying a rescaling of the loop momenta k ^ Xk does not change the value of the 
integral, and we thus have 

Jd'k^^^O. (3.13) 
The same holds for integrals of the form 

since by dimensional analysis (the result must have mass dimension d — 2 {rii + ^2), 
and p1 is the only Lorentz invariant quantity out of which a mass scale can be 
formed) they must be proportional to 

an expression that vanishes for Ked > 2 (rii + 722). Thus, via analytic continuation, 
they vanish in the entire complex d-plane. Additional scalar products in the numer- 
ator do not change this, so we may identify entire topologies (diagrams with a given 
set of propagators) as vanishing. The generahsation to three-loop integrals brings 

^This statement is true for any two- or higher-loop (one-particle-irreducible) vertex graph, 
definition of tadpole and snail topologies is given in the Qgraf manual, see [4]. Loosely 
speaking, a tadpole contains a propagator that carries zero momentum. Fig. 3.5 contains a snail 
that is not a tadpole. 
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about a complication, namely that the three-loop integral might factorise into a 
1-loop and a 2-loop integral (or three 1-loop integrals). In that case, each factor 
must be checked separately. 

Building upon the results obtained in the previous section has one main advan- 
tage: It is easy to factorise the integral, because in our given basis, the number of 
propagators that contain more than one loop momentum is minimal. The integral 
thus factorises if out of the propagator momenta ki — k2, ki — k^, k2 — k^ at most 
one is present in the given topology. Furthermore, since in our basis no propaga- 
tor depends on p2 alone, at least one propagator has to be of the form ki — q and 
thus introduce the only nonvanishing scale into the integral. If an integral only 
contains factors ki — q and kj — pi, a. shift in the loop momenta ki ki+pi shows 
that the integral only depends on p2 and thus vanishes. To summarise, the following 
algorithm detected most vanishing planar topologies: 

1. Bring the integral into a standard form, with propagator momenta out of the 
set of the auxiliary topology. 

2. Factorise the 3- loop integral into lower- loop integrals as far as possible. 

3. For a topology to be non- vanishing, require each factor to contain at least one 
propagator 1/ {k^ — qf and one propagator 

One type of topology that vanishes but was not detected as such using the above 
algorithm is shown in fig 3.5. This particular diagram is harmless since Qgraf has a 
switch to turn off snail or tadpole diagrams. On the other hand, this example shows 
that shifts in loop momenta play an important role in the above algorithm, since 
after a shift in the loop momenta, this topology 

{ki, ki - pi, ki -q,ki- k2, ki - /cg, ki - ks} 

k2 -A';2 + A-i, 
ks — >■ —ks + ki 

> {ki, ki - pi, ki - q, k2, k^, - {ki - ks)} (3.15) 

will be identified as vanishing since the factor f (i^/corf^fc^ ,2/,'^ — rr^ has no scale. 
This illustrates again that the propagators that involve two different loop momenta, 
i.e. ki — k2, ki — ks, k2 — k^ should be eliminated if possible such that the factorisation 
of the integral in the above algorithm is maximal. Which shifts should be tried 
is not clear (we will discuss this point in the next section). Here, we state that 
even without any further shifts, all vanishing topologies except the one mentioned 
above were detected. As we will describe in sec. 4.3 one may double-check using 
integration-by-parts identities. The two methods yield identical results for diagrams 
that do not contain snails. 
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Figure 3.5: A vanishing planar diagram. It contains the propagators ki — k2, ki — k^, 
/c2 — k^ and e.g. ki — q. Looking at the propagators, it is easily seen to be vanishing 
only after a shift k2 —>■ k2 + ki, k^ —>■ ks + ki is applied. 

3.4 Exploiting Symmetries between Topologies 

If two Feynman integrals are equal, there is usually an easy way to understand this 
in terms of a symmetry of the integrand. These symmetries could be: 

• a permutation of the labels of the loop momenta ki k2 ^ k^, 

• the interchange of pi <-> P2, 

• a shift of the loop momenta. 

It should be clear that the first two operations actually transform entire topologies 
into other topologies (the explicit exponents of the propagators do not matter for 
this transformation). Shifting loop momenta is more delicate since it is not clear 
which shifts of loop momenta might be useful, and because not all shifts transform 
entire topologies in the same way. For this reason we discuss shift and permutation 
symmetries separately in the next two sections. 

3.4.1 Symmetries that are Permutations of Momenta 

Permutations of loop momenta correspond to a mere relabelling of the integration 
variables, i.e. the loop momenta. Interchanging pi and p2 cannot alter the value 
of (massless) vertex integrals since = 2 pi ■ p2 is the only scale left after the 
integration has been carried out. In the planar auxiliary topology, the permutation 
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Pi ^ P2 can be performed by applying the shift ki — > —ki + g for all i = 1, 2, 3 
simultaneously, followed by the permutation pi <-> p2- It implies 

Pi ^ P2 : 

T(ni, n2, ns, n^, n^, uq, n^, ns, rig, nio, nn, ni2) = (3.16) 
T(ns, nio, n^, n^, n^, Uq, n?, Ui, Ug, n2, ^11,713) . 

Similarly, we have 
/ci <-> A;2 : 

T(ni, n2, ns, n^, ^5, rie, 717, ng, rig, riio, nn, ^12) = 

T(n2, ni, na, 77,4, ne, ^5, ng, nio, ny, ng, nn, 77,12) , 

ki ks : 

T(ni, 77,2, na, 714, 715, uq, nr, ng, ng, nio, nn, ni2) = (3.17) 
Tins, n2, ni, ng, n^, rii, nn, ni2, ng, nio, ny, ng) , 

k2'r^ k'i-. 

T(ni,n2, ns, n4, n^, Uq, nr, ng, ng, nio, nn, ni2) = 

T(ni, ns, n2, ns, n4, ne, ny, ng, nn, ni2, ng, nio) • 

Clearly, such relations should be determined and topologies related by them should 
be identified, reducing the number of topologies that have to be evaluated. 

If two topologies have such a relation, we need a prescription which of the two 
topologies should be mapped onto the other. An ordering of topologies is provided by 
their topology ID number, given as the number & {n ^J^,n < 2^^ — 1} whose binary 
representation simply indicates which of the 12 propagators out of the set (3.5) are 
present in this topology. All integrals in this topology then receive the same ID 
number. E.g. for the topology (3.7) (see fig. 3.1) we write 

T(-1,-1,0,-1,-1,-1,-1,0,0. -L-1,-1) 

^ ID = 1 + 2 + 8 + 16 + 32 + 64 + 512 + 1024 + 2048 = 3707 = 1110011110116 . 

(3.18) 

We will add a corresponding subscript to the planar topology T = TJd- Imple- 
menting permutations of loop momenta is now straightforward. In this example, 
permuting ki k^ yields 

TsM-i, -1, 0, -1, -1, -1, -1, 0, 0, -1, -1, -1) 

ri7go(0, -1, -1, -1, -1, -1, -1, -1, 0, -1, -1, 0) (3.19) 

with an ID = 1790 < 3707. The permutation pi ^ p2 corresponds to the shift 
ki > ki + q,pi ^ p2 for aU i = 1, 2, 3 simultaneously, which in this case does not 
lower the ID number. 
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3.4.2 Symmetries that are not Permutations 

There are also more complicated permutations that one might consider trying. For 
example the topology with ID = 187 has the propagator momenta 

ki, k2, ki - k2, ki - ks, k2 - k^, ki - q . 

Applying the shift 

fci ^ k[ 

k2 ^ + ki (3.20) 
^3 -k'3 + ki , 

one ends up with topology 175 that has momenta 

ki, k2, k3, ki - k2, k2 - k^, ki - q . 

But under the shift (3.20) additional scalar products such as k2—pi, k^—pi transform 
into a linear combination of inverse standard propagators, and thus integrals in the 
topology 187 are expressed as several terms in topology 175. Since this might enlarge 
the size of integration-by-parts identities (see sec. 4.1), it is not clear whether they 
are of any use there. We did not use these identities at all, although they might be 
of moderate use to reduce the number of subtopologies that appear as integrals in 
the form factor. 

We have not answered the general question of what groTip of (linear) loop mo- 
mentum shifts would map subsets of the auxiliary topology 1 (the propagators of a 
topology) into subsets of the same set (3.5). If the sets of propagators are different, 
it seems that only in rare cases there exists a shift that is not a permutation of 
loop or external momenta. One of the few examples was given in eq. (3.20), another 
in fig. 3.5. The far more frequent case is that a shift permutes the propagators 
of a given topology. This happens for example whenever the integral contains a 
two-point function (also called a bubble insertion) as in fig 3.6. There, the shift 
ks — >• —ks + ki interchanges the two indicated propagators and all other propagators 
are left unchanged (but not the scalar products). Thus we obtain what is called a 
symmetry identity 

T734 - ^734|fc3^_fc3+fc^ = (3.21) 

The significance of symmetry identities will be discussed in sec. 4.4. Basically, if a 
symmetry is present, it is possible to gain a lot of relations between integrals with 
little computational effort. Since these symmetry identities are found by inspec- 
tion of the diagram anyhow, it was important that the symmetry considerations 
in sec. 3.4.1 together with the algorithm in sec. 3.3 were enough to bring down 
the number of different topologies to about 200. Their symmetry properties were 
then implemented by hand. Drawing those diagrams, no relations except the ones 
mentioned here were apparent. 
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Figure 3.6: The topology with ID=734 has a bubble insertion and thus an additional 
symmetry: The two propagators with momenta and ki — k^ can be interchanged 
with a shift ks — > —ks + ki. 



Chapter 4 

Methods for the Reduction 



In the previous chapters we have generated the three-loop QCD quark form factor 
T {q"^) as a Unear combination of three-loop topologies T. Obvious symmetries 
have been used and all vanishing topologies were removed. At that stage, our 
representation of the form factor consisted of 2951 integrals. Our next main task 
will be to use integration-by-parts identities in an automated way to relate these 
integrals to a small number of master integrals. The computation of some of these 
master integrals is very involved and not part of the present work. By now most of 
the master integrals have been computed [28, 29] , such that the value of the form 
factor (it is just a number) will be in reach once the reduction of the form factor to 
master integrals is achieved. 

4.1 Integration-By-Parts and Lorentz Invariance 
Identities 

Integration-by-parts identities [30] (IBPs) are relations that hold generally in dimen- 
sional regularisation. They express the fact that dimensionally regulated integrals 
are Poincare covariant. Consider a general integral / depending on loop momenta 
ki . . .kn' For an infinitesimal translation ki ^ ki + vj, we have 



/ [ki . . . , ki, . . . , kn) 



I{ki 



. , ki -\- Vj, . . . , kn) ■ 



(4.1) 



Expanding in Vj, we get 



^ dkf 



I{ki 



(4.2) 



Integrating by parts the IBP identities in our case read 




(4.3) 
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with Vj e {ki,k2,k3,pi,q}, and where I denotes the integrand of a three-loop (mass- 
less vertex) integral /. 

Acting with the derivative in (4.3) on ?;^X generates a number of terms, with 
one exponent of a propagator incremented by one and with contracted with a 
four- momentum. The scalar products can be rewritten using (3.6), yielding a linear 
combination of functions T. In this way IBP identities relate T's with different 
(integer) arguments. Thus IBP identities form a system of equations that are linear 
in the subtopologies T. 

We call T a seed integral since it generates the IBP identities. Below, we will often 
call the integrals that appear in the IBP identities of a seed integral its subtopologies. 
Note that -^k^ = ci in dimensional regularisation, so that the coefficients in the 
system of IBP identities contain the symbol d — A — 2e. 

It has soon been realized that the number of IBP identities grows faster than 
the number of unknowns. This apparently overconstrains the system. In fact the 
system is underconstrained, so that its solution relates the subtopologies it contains 
to so called master integrals. These are integrals with a minimal sum of exponents 
of propagators r and sum of exponents of scalar products s that cannot be reduced 
to integrals with a smaller number of propagators t using IBP identities. In our 
case, we expect all master integrals to have r = t, s = 0. In section 4.2 we describe 
an algorithm [31] to use IBP identities in an automated way to express more com- 
plicated integrals as a linear combination of simpler ones, and eventually in terms 
of master integrals. 

Just as IBP identities express the translation invariance of dimensionally reg- 
ulated integrals, we may also exploit Lorentz invariance in the same manner [32]: 
An infinitesimal Lorentz transformation on an external momentum reads — >■ 
Pi + ^Pi = Pi ^^uPi with Se'^ = —Se'^. Again, expanding the identity 

I (pi, . . . ,Pn) ^ I {Pl + Spi, . . . ,Pn + Spn) (4.4) 

we obtain 

or, using the antisymmetry of 5e'^, 

For three-point functions, we have n — 2, i.e. 

By contracting these equations with antisymmetric combinations of pi, . . . ,p„ we 
get the Lorentz invariance (LI) identities. For n — 2, the only such combination is 
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- P2,M and we get 
or with q = pi + P2, 

LI identities have recently been proven to be hnearly dependent on IBP identi- 
ties [33]. While it is possible to obtain these relations from IBP identities, this is 
computationally expensive because a larger system of IBP identities must be solved. 
For that reason we used LI identities at all stages of our reduction. 

The question of which integrals have to be used to generate IBP and LI identities 
will be discussed in sec. 4.4. 

4.2 The Laporta- Algorithm 

In [31] an algorithm (by now simply called the Laporta algorithm) was proposed 
that makes use of IBP (and LI) identities to reduce integrals to master integrals in 
an automated way. The algorithm solves a system of IBP identities using a specific 
ordering of the integrals according to their complexity. To quantify this, we now 
define the three numbers mentioned in the previous section that characterise an 
integral T (ni, . . . , ni2) (n^ G Z): 

t: the number of propagators: t — J^^^i {—n-i — 1) 

r: the sum of the exponents of the propagators: r = ^1=1 {i^-il (— rij) 

s: the sum of the exponents of the scalar products: s = \ni\6 (ui) 

(where ^ (x) = 1 for a; > and otherwise; remember that propagators have < 0) . 
The Laporta algorithm instructs us to order integrals first according to their t, then 
their r, and finally their s. Integrals with the same t, r, s will be further ordered 
according to the number of terms in their coefficient, and then according to the total 
number of terms in their equation. 

The system of IBPs may then be solved using Gauss' algorithm (further detail 
on the algorithm are given in sec. 4.6). All the identities extracted from the system 
of IBP identities should map more complex integrals (according to the ordering just 
given) to simpler ones. If enough identities are solved, the solutions should express 
all desired integrals as linear combinations of master integrals only. 
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4.3 Identifying Subtopologies 

It is a good idea to look ahead and try to estimate the scale of the calculation we 
are aiming at. The form factor has already been introduced and expressed as a 
linear combination of three-loop integrals. We will next try to estimate the size of 
the system of IBP identities that has to be solved in order to reduce the integrals 
in the form factor to master integrals. The master integrals themselves will also be 
identified in the process. 

To start we should identify the topologies that might arise when we reduce 
our form factor. Here, we do not yet care about which integrals appear in the 
reduction. We will have to perform a reduction for the integrals in these topologies 
by solving IBP identities. Of course, equivalent topologies (e.g. those that differ by 
a permutation of loop momenta) only have to be reduced once since we know the 
equivalence classes of topologies (see sec. 3.4 and the appendix B). 

The t — 1 subtopologies for a topology with t propagators arise in the following 
way: In an IBP identity a scalar product is formed by the contraction of with a 
vector in the numerator of When it is expressed as a linear combination of 

inverse propagators, an inverse propagator may cancel one of the existing propaga- 
tors. Thus subtopologies arise by removing one propagator from the given topology^. 
Which of the propagators are cancelled can only be determined by writing down the 
IBP identities. However, we need to write down the IBP identities for the topology 
with r — t and s = only: Any possible cancellation of a scalar product against a 
propagator will show up here, and therefore alH — 1 subtopology ID numbers can 
be determined in this way. 

Two interesting properties of the subtopologies can easily be derived from the 
IBP identities of the r — t, s — seed integral of a topology: Firstly if the topology 
is vanishing, this will almost always be apparent after the system of the 15 IBP 
identities from the r = t, s = seed integral is solved. Secondly, if the topology is 
reducible, the reduction equation for the r = t, s = integral can often (but not 
always) be extracted from the same solution of 15 IBP identities. Thus we arrive 
at a natural procedure to get a first insight into the subtopologies that have to be 
reduced: 

1. Create a list Ljd of the ID numbers of all topologies present in the form factor. 

2. Set ^0 equal to the maximum number of propagators present in the topologies 
in Lid; in our case, to — 9. 

3. For every topology in Ljd with t — to, generate the IBP identities for its r — t, 
s — seed integral and 

• store all subtopology ID numbers found in the identities and add them 
to Lid; 



In symmetry identities, up to s propagators may be cancelled. 
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f 

t 


■fj- \J1. ±J_/o 


•jf- U± 11 1 CU.. ±J— 'O 


S of "hard" IDs 


9 


3 


1 


1 


8 


25 








7 


46 


1 


4 


6 


33 


3 


7 


5 


10 


4 





4 


1 


1 






Table 4.1: The number of topologies and irreducible topologies found as subtopolo- 
gies in the three-loop quark form factor. The third column shows the number of 
topologies that have to be reduced numerically. Symmetries are already taken into 
account. 

• solve the system, and determine from the solution whether the topology 
is vanishing or else whether it is reducible or (possibly) irreducible. 

4. Set to = to — 1 and repeat the previous step. 

This creates a hst of all (sub-)subtopologies that will occur during the reduction. 
When generating the IBP identities, the information on the equivalence of certain 
topologies derived in chapter 2 should be used. 

One subtlety should still be pointed out: Several topologies in fact are reducible 
although they are not found to be so using the above algorithm. One example is 
the topology with ID=734 (fig. 3.6). Their reduction equations are found only after 
the inclusion of additional IBP identities, e.g. uptor = t-|-l,s = l. Since these 
systems are already quite big, it makes sense to check reducibility first using the 
r = t, s = seed integral only, and checking the potentially irreducible graphs in a 
second run using e.g. t<r<t+l, s<l seed integrals. 

Table 4.1 gives the number of nonvanishing subtopologies and irreducible subtopolo- 
gies found in the appfication of the above algorithm and after double-checking the 
irreducibility of the topologies. Obvious symmetries (such as permutations of the 
loop momenta, see sec. 3.4) have already been taken into account. Out of the topolo- 
gies, some can be integrated loop-by-loop or reduced using symbolic IBP identities 
and do not pose a problem (see chap. 5 and sec. 4.5). All others (called "hard" here) 
are shown in fig. 4.1 ior t — 6 and in fig. 4.2 for t — 7. They have to be reduced 
using explicit numeric IBP identities^. The irreducible master integrals for t — 4,5 
arc given later in fig. 5.2. 

•^Notc though that the first four topologies in fig. 4.1 {t = 6) and the first three topologies in 
fig. 4.2 (t = 7) contain a two-point function inserted into the same two-loop topology. One option 
that we have not yet fully explored is to reduce the two-loop topology with one exponent of a 
propagator kept generic. This introduces the symbolic exponent as a second symbol (apart from 
the space-time dimension d), but reduces the number of equations per seed integral from 16 to 7, 
and also the number of seed integrals per r-s-block is reduced considerably. 
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ID=691 

Figure 4.1: The "hard" 6-propagator topologies appearing in the reduction, i.e. 
those that are not integrable loop-by-loop or for which no symbolic reduction equa- 
tions were found. The topologies are denoted using their binary ID number as well 
as (if present) the name used in [28]. They spht into two groups: There are the 
irreducible topologies (the middle column, consisting of the topologies 429, 444 and 
691) that are not integrable loop-by-loop. The second group (first and third col- 
umn) consist of reducible topologies for which no symbolic reduction equations were 
found. All these graphs have to be reduced using explicit numeric IBP-identities. 
Note that the first four and also the last two topologies are bubble-insertions into 
the same 2-loop topology. 
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ID=758,A7,3 



Figure 4.2: The "hard" 7-propagator topologies appearing in the reduction. The 
topologies are denoted using their binary ID number as well as (if present) the name 
used in [28]. They should be reduced using IBP-identities where the exponents of 
propagators are written out explicitly. Note that the topologies in the first row are 
bubble-insertions into the same 2-loop topology. They are all reducible. 
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The result that there is only one topology left ai t — A is expected: Generally, 
for any diagram the relation 

L-l^t-V (4.10) 

holds, where L is the number of loops and V the number of vertices. Diagrams with 
only one vertex are snail diagrams and vanish (in massless, dimensionally regulated 
computations). Thus t > A follows. Moreover, for t = 5, L = 3 it follows that V = 3. 
But V = 3 means that all integrals contain at most one three-point function. They 
can be integrated loop-by-loop, as will be discussed in chap. 5. For t > 6 the 
integrals become complicated, and we have to resort to the Laporta algorithm to 
reduce integrals to master integrals. 



4.4 A First Estimate of the Scale of the Compu- 
tation 

We have now gained an overview of which topologies have to be reduced using the 
Laporta algorithm. We next ask which t-r-s-classes of integrals might appear in our 
reduction. 

Given a seed integral It,r,s for the massless quark form factor, the derivative 
acting on a denominator will increment r by one. At the same time, the vector 
y/j- from the IBP identity will be contracted with a vector that was obtained from 
the differentiation via the chain rule. This yields an additional scalar product in 
the numerator. The results are integrals of the type^ Tt^r+i,s+i or Xj^^+i,* • When 
expressed as a linear combination of inverse propagators, the scalar product may 
cancel with one propagator, yielding either g or Xt_i,r,s integrals. Finally, integrals 
of the type 1t,r,s-i may be generated when the derivative acts on a scalar product 
in the numerator. It may remove one loop momentum from the scalar product and 
substitute it by an external momentum. To summarise. 



generate IBPs 

t, r, s > < 



' t,r + l,s + l 
t, r + l,s 

t,r,s . (4.11) 

t,r,s — 1 
t — l,r,s 



The topologies 2t-i,r,s are the simplest in this list and will appear in the final re- 
duction equations. The important observation here is that if a seed integral Tt^r,s is 
used to generate IBP identities, this introduces subtopologies {t — 1 integrals) with 

the same r and s. 

Concerning possible subtopologies, symmetry identities are quite different from 
IBP identities. They mostly relate integrals that have the same t, r and s as the 



^In case the scalar product is pi .p2 we do not count this as a scalar product that increments s 
since it factors out of the integral and does not increase the complexity of the integral. 
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seed integral. This is a significant advantage: To obtain relations for integrals with 
a given r and s, IBP identities are generated using mostly the seed integrals with 
r — 1, s — 1, whereas the symmetry identities are generated using the seed integrals 
with the same r and s. Clearly, there are many more seed integrals of the type r, 
s than of the type r — 1, s — 1 orr — 1, s. Symmetry identities thus provide a 
significant amount of information, even though there is mostly a single symmetry 
identity present, compared to the 15 IBP identities. Indeed we will see shortly that 
there is roughly an order of magnitude more r + 1, s + 1 than r, s seed integrals, 
such that if present, a symmetry identity provides a number of relations between 
r-s integrals comparable to the one of IBP identities. 

After the generation and basic manipulation of the three-loop quark form factor, 
we find by inspection that the most difficult integrals to be computed have 



The corresponding topologies finally reduce to t > 6 topologies that are irreducible, 
or to t = 5 topologies that can be integrated loop-by-loop. 

Next, we compute the number of seed integrals of a given topology as a function 
of t, r and s. This is simple combinatorics, with the result [32] 



Since we have 15 IBP identities and one LI identity, the corresponding number of 
equations for seed integrals with a given t, r and s is just Nt^r,s — ^^^^ {^t,r,s)- The 
explicit numbers relevant for us arc given in table 4.2. Since we will evaluate the 
integrals explicitly for t = 5, table 4.3 gives the number of seed integrals in a given 
r-s-block, Nri^r,s- 

The numbers are not very pleasing. The code Solve that will be discussed in 
sec. 4.6 will turn out to be able to solve systems with up to ~ 3000 equations. In 
sections 4.7 and 4.8 we will improve this code in various ways. This will allow us 
to solve systems with close to 10000 equations. But clearly, we cannot afford to 
reduce the topologies with t, r, s as in eq. (4.12) using seed integrals up to r = 9, 
s = 4 for f = 7, 8, 9. We would have to expect subtopologies with t — 6, r — 9, 
s — A. To reduce those topologies, we expect it to be necessary to solve a system of 
equations with an order of magnitude more equations than could be handled using 
Solve: Tab. 4.2 predicts 188160 equations to be solved in the r = 9, s = 4 block 
alone. The fact that thesis does not present the full reduction of the three-loop 
quark form factor shows that ultimately we have not been able to fully address this 
challenging problem. We will discuss various aspects of it including some attempts 
to solve such large systems in the remainder of this thesis. 

On the other hand, for 9- (and 8-)propagator topologies the number of equa- 
tions is quite moderate. We may go ahead and reduce these systems without much 
difficulty, given that we have a CPU cluster at our disposal. 



t^9 
r = 9 
s = 4 



(4.12) 




(4.13) 
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t = 9 




s=0 


s=l 


s=2 


s=3 


s=4 


r=9 


16 


61 


160 


320 


560 



t = 8 




s=0 


s=l 


s=2 


s=3 


s=4 


r=8 
r=9 


16 
128 


80 
640 


240 
1920 


560 
4480 


1120 
8960 



t = 7 




s=0 


s=l 


s=2 


s=3 


s=4 


r=7 


16 


96 


336 


896 


2016 


r=8 


112 


672 


2352 


6272 


14112 


r=9 


448 


2688 


9408 


25088 


56448 



t = 6 




s=0 


s=l 


s=2 


s=3 


s=4 


r=6 


16 


112 


448 


1344 


3360 


r=7 


96 


672 


2688 


8064 


20160 


r=8 


336 


2352 


9408 


28224 


70560 


r=9 


896 


6272 


25088 


75264 


188160 



Table 4.2: The number of equations generated from a given r-s-block of seed integrals 
for a system with 12 independent scalar products. It grows quickly with r and s 
and with decreasing t. 







t 


= 5 








s=0 


s=l 


s=2 


s=3 


s=4 


r=5 


1 


8 


36 


120 


330 


r=6 


5 


40 


180 


600 


1650 


r=7 


15 


120 


540 


1800 


4950 


r=8 


35 


280 


1260 


4200 


11550 


r=9 


70 


560 


2520 


8400 


23100 



Table 4.3: The number of seed integrals in a given r-s-block for a system with 12 
independent scalar products. 
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4.4.1 The Hierarchy of Integration- by-Parts Identities 

The most immediate objection to the above argument that the system of IBP iden- 
tities will grow unmanageably large is that to reduce an integral It,r,s, there might 
be no need to generate and solve IBP identities from seed integrals with the same 
values for t, r, s. As is apparent from (4.11), a seed integral Itjr~i,s-i or Xt^r-i,s 
also generates identities involving It,r,s, and their t — 1 subtopologies are of the form 
It-i,r-i,s-i and Xt_i s, respectively. This is illustrated again in fig 4.3. For the 
reduction of the integrals in (4.12), t = 9, r = 9, s = 4 seed integrals will be neces- 
sary. This will lead to subtopologies with t — S, r — 9, s — A. For their reduction, 
one may use the following seed integrals: 



seed integral 


r=9.s= 1 generated in 


t = 7 subtopologies 


t = 9, r = 9, s = 4 
t = 8, r = 8, s = 3 
t = 8, r = 8, s = 4 
t = 8, r = 9, s = 4 


t — 1, r, s 
t, r + 1, s + 1 
t, r + 1, s 
t, r, s 


none 

r = 8, s = 3 
r = 8, s = 4 
r = 9, s = 4 



The first line corresponds to so-called left-over equations. These solutions of the 
system of IBP identities contain only integrals of type t—1, r, s. They might 
relate different topologies and reduce the number of t = 8, r = 9, s = 4 topologies 
needed to be reduced by other identities. The second line corresponds the dashed 
blue arrow in fig. 4.3. These identities are more useful than the ones in the third 
line (corresponding to the solid black arrow in fig. 4.3) since they do not contain 
r = 8, s = 4 7-propagator subtopologies. Finally, identities using t = 8, r = 9, s = 4 
seed integrals should be avoided in any case since they lead tot = 7, r = 9,s = 4 
subtopologies. 

To summarise, in our context it is crucial to note that some IBP identities are 
more useful than others, since they contain less complicated subtopologies. If several 
identities reduce the same integral, one must use the identity that was generated 
from seed integrals with the lowest possible r and s. In practice, this means that 
different collections of the same set of IBP identities must be solved and substituted 
in. For example one should use the solution for the system with t — 8, r — S, 
s < 3 before using the solution for the system with t = 8, r = 8, s < 4. A routine 
dedicated to finding in such an incremental way the solution of a system of IBP 
identities will be presented in sec. 4.8.2. 

4.5 Using Symbolic Integration-by-Parts Identi- 
ties 

As a result of the above investigations, we find that all r = 9 and all s = 4 integrals 
could be reduced to t — S, r — S, s — A or t — 8, r = 9, s — 3 integrals. However, 
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t=8 



8 








Figure 4.3: Some integration-by-parts identities are more useful than others: To 
reduce integrals with t = 8, r = 9, s = 4, we may use seed integrals with e.g. r = 8, 
s = 3 (dashed blue arrow) or r = 8, s = 4 (solid black arrow). But r = 8, s = 3 
seed integrals lead to IBP identities without any t = 7,r = 8, 3 = 4 subtopologies, 
whereas IBP identities from r = 8, s = 4 seed integrals do contain these more 
complicated subtopologies. 



using IBP identities from seed integrals with t = 8, r = 8 and s < 4 in the way 
described in the previous section, it was not possible to reduce all 8-propagator 
integrals. 

As a next step, one could include t = 8, r = 9 seed integrals into the system of 
IBPs, producing t = 7, r = 9 subtopologies. This consequence cannot be avoided 
even after the inclusion of symmetry identities. However, as we will show now, the 
fact that all 8-propagator topologies are reducible (see tab. 4.1) leads to an elegant 
way of reducing them using just 15 symbolic IBP identities and one LI identity (per 
topology). The subtopologies will even tend to be simpler when using symbolic IBP 
identities as compared to using "numeric" IBP identities with higher r and s. 

By symbolic IBP identities, we mean IBP identities that are written down for 
general exponents ni,...,ni2. Symbolic integration-by-parts identities have been 
used already in [30]. It is perhaps easiest to illustrate their use in the case of the 
one-loop vertex-correction, fig. 1.1. The three-propagator integrals that appear in 
its evaluation are of the form 



2^1-ioop (ni,n2,n3) = J d'^k- 



[kr- {{k - p,fr {{k - q)T ' 

Again with q = pi + P2- When writing down the IBP identities for this case, 

d 



dki" 
d 

d 
dki" 



FJi_loop (?^l,?^2,^3) = 



(4.14) 

(4.15) 
(4.16) 
(4.17) 
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it is useful to introduce the following operators: 

l+Ji_ioop (ni, 712, ns) = Ji_ioop {ni + 1, n2, ris) , (4.18) 
2+Xi_ioop {ni, 712, ris) = Xi_ioop (ni, n2 + 1, ng) , (4.19) 
3+Xi_ioop (ni, 7i2, na) Xi_ioop {rii, 7i2, ria + 1) , (4.20) 

and similarly for 1~, 2~, 3~. Using this notation the three IBPs (4.15), (4.16), and 
(4.17) read 



(7l2l"2+ + 7l3l"3+ - 7l3g^3+ 

- d + 2ni + n2 + rzs) Ji_ioop {ni,n2, ns) = 
(-n2l"2+ - n3l"3+ + nsg^S^ 

+n32~3^ + nil^2~ -ni + n2)li-ioop{ni,n2,n3) = 
(rzil+2~-nil+3" + nig^l+ 

-ri22+3" + ri32"3+ + ri2 - najXi.ioop (rii, 7i2, ^3) = 



(4.21) 
(4.22) 
(4.23) 



The most complicated integrals in (4.21) and (4.22) are contained in those terms 
that do not involve a lowering operator i~, i.e. the terms ^q^n^S'^. We eliminate 
these terms by adding (4.21) and (4.22). We get 

(7132-3+ + 7iil+2- - d + 7ii + 2712 + ns) Ti_ioop (ni, 712, ria) = . (4.24) 

For general rii, all terms seem equally complicated, and the Laporta algorithm would 
thus instruct us to solve for the first or second term, since their coefficients are small. 
But this would introduce inverse coefficients 1/ni or 1/713 which would actually 
become singular when the reduction of the topology should take place. Indeed, 
solving for the second term in (4.24) and relabelling the Tij's, one gets 

T ( . d-ni-2n2-n3-l , . 

2-l-loop (ni, n2, n-i) = Il-loop (^1 - 1, 712 + 1, ^3) 

7Zi - 1 (4-25) 
- n-i Xi-ioop (7ii - 1, 7i2, na + 1) 

which cannot be evaluated for n\~\. 

Instead, eq. (4.24) should simply be solved for Xi_ioop n2, n^). We may write 
this as an operator equation 

(711 1+2- +7132-3+) , (4.26) 



d — rii — 2n2 — 
or more explicitly as 



2^1-ioop {ni, 7i2, ns) = -z ^- (711 Xi_ioop {ni + 1, 7i2 - 1, ria) 

d - 7ii - 27i2 - ria \ (4.27) 

+ 7l3Xi_ioop (7li, 7l2 - 1, ^3 + 1)^ . 
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Here, the fact that the inverted coefficient contains the symbol d prevents the denom- 
inator from becoming singular. For irreducible topologies, the inverted coefficients 
do not contain the symbol d and thus the equation is not applicable to the master 
integral. 

Even for reducible topologies, it is by no means obvious how symbolic identities 
may be used in general. In the above example it was possible to combine two 
equations to eliminate a term that does not have an operator i~, n^q'^S'^. This 
is not always possible, and so some reducible topologies have to be reduced using 
numeric IBP identities. 

To generalise this to three-loop (reducible) integrals, IBP identities arc solved 
just as in the Laporta algorithm^ until topologies with the same r and s as the seed 
integral are reached. At that point, integrals with coefficients containing a symbol 
d should be favoured. The corresponding integral usually has no lowering or raising 
operators acting on I. 

The symbolic identity found in this way may then be applied repeatedly until one 
propagator is cancelled completely, and the topology thus changes. This substitution 
has to be detected and a further apphcation of the symbolic identity has to be 
prevented. The next section will report in more detail on the reformatting necessary 
for symbolic IBP identities. 

4.5.1 Symbolic Integration- by-Parts Identities for Irreducible 
Topologies 

So far the experience gained with symbohc IBP identities for reducible topologies 
has not let to any significant insight in how symbolic IBP identities could be used for 
irreducible topologies. Symbolic IBP identities for irreducible topologies always seem 
to have coefficients that would become singular when applied to certain integrals, 
especially to the master integral. When they are applied recursively to a general 
integral, some exponents might be lowered until ultimately they are not applicable 
any more since that would again lead to singular coefficients just as in (4.25). 

A second caveat lies in scalar products (i.e. propagators with non-positive expo- 
nents) whose exponents are decremented upon application of the symbolic identity. 
If no scalar product is present in the integral onto which the identity is applied, an 
additional propagator is generated, which is clearly not desirable. So the application 
of the identity to those cases should also be prevented, except if the coefficient of 
those integrals becomes zero. In this case the identity may still be applied. 

Still, for sufficiently large exponents rii some of these identities do contain valu- 
able information in that they are indeed able to reduce r and s for certain integrals. 
Generating symbolic IBP identities is computationally cheap, and applying them is 
worth trying at least in those cases where otherwise large r-s blocks of IBP identities 
would have to be evaluated. But most of the identities actually are not useful at 



For the evaluation of r and s, set all Uj equal to unity. Drop the ordering according to t. 
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all, in that they relate a given integral of a certain complexity to many integrals of 
the same complexity. Thus symbolic identities should only be used if the identity, 
when applied, in fact reduces r or s of the integral it was applied to. 

From the symbolic identities contained in the solution of the system of (symbolic) 
IBP identities, the ones must be selected that match the above criteria. All other 
identities arc discarded. The identity then has to be rewritten as FORM code in 
such a way that it is only applied when appropriate. Altogether our code that 
analyses the identities is rather complex. It was thus written in the more versatile 
language of Mathematica [34]. 

The following sample code illustrates the final form of the identities (we have 
introduced ellipses to condense the code): 

1 id Topol(-7,1694,[] ,ml?,n2? ,ml2?) = 

2 deltap_(l + n2)*Topol( -7 ,1694 ,[] ,ml,n2 , . . . ,ml2) 

3 + delta_(l + n2)*Topo2( -7 ,1694 ,[] ,ml,n2 , . . . ,ml2) ; 

4 id Topol( - 7 ,1694 ,[] ,ml?posO_ , . . . , n8?neg0_ ,0 ,nlO?negO_ ,...) = 

5 Topo2( -7,1694 ,[], ml, .. . , n8 , , nlO , . . . ) ; 

6 id Topol(-7,1694,[] ,ml?,n2? ,ml2?) = 

7 - ml/(l+n2)*Topol(-7,1694,[] ,-l+ml,H-n2,... , -l + m9, . . . ,ml2) 

8 - n3/(l+n2)*Topol( -7,1694,[] ,ml,l + n2 ,. . . ,ml2) 

9 + Topol( -7,1694 ,[], ml, l + n2 ,ml2) 

10 ... 

11 + ml2/(l + n2)*Topol(-7,1694,[] ,ml,l + n2,... ,-l+ml2) 

12 + nlO/(l + n2)*Topol ( -7 ,1694 ,[] ,ml,l + n2 , . . . ,ml2) 

13 + n8/(l+n2)*Topol( -7 ,1694 ,[] ,ml,l + n2 , . . . ,ml2) ; 

14 id Topol(-7,1694,[] ,ml?,0,n3? ,ml2?) = 

15 Topo4 ( - 7 , , [] ,ml , , n3 , . . . , ml2 ) ; 

The denominator of this identity would become singular if evaluated at ^2 = — 1. 
Thus lines 1-3 prevent the identity to be evaluated if 712 = —1 by renaming Topol 
to Topo2 if n2 = —1. 

On the right-hand-side of the identity (line 7) the scalar products mi, nig and 
mi2 are lowered, but only when rrig is lowered from to -1 there is a danger of 
introducing a higher-i topology, since the coefficient of the other integrals contain 
a factor mi and mi2, respectively. They make these terms vanish if mi = or 
mi2 = 0. Thus two lines are written (lines 4-5) that prevent the case where mg is 
lowered, but no such code is written to prevent the lowering of mi or mi2. 

Finally, lines 14-15 rename topologies where one propagator disappeared, as is 
the case with the exponent of the second propagator (— ^2) that is lowered in this 
identity, 712 — > ^2 + 1. In this way the whole code block can simply be iterated a 
few times until it is not applicable any more. 

Also, note that in fact this identity has correctly been selected as useful: all the 
right-hand-side topologies have the exponent of the second propagator (—1 — 
lowered by one with respect to the left-hand-side, a property that the symbolic IBP 
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identities (right after their generation) usually do not have. It only arises in some 
of the equations that constitute the solution of the system of IBP identities. 

The identity as it stands in this form can be safely applied to the form factor 
expression that should be reduced. With some luck, high-r-s subtopologies can 
be reduced using these identities. But mostly integrals do not find a reduction 
equation because in our appUcation (t < 9, 6 < r < 9, s < 4) most exponents of 
propagators and scalar products are still 1 and 0, respectively, and because symbolic 
IBP identities for irreducible topologies TisTially become singular at some rii = — 1. 
Therefore it is quite likely that for most integrals none of the identities matches all 
criteria listed above. 



Solve is the name of a code by E. Remiddi that implements the Laporta algorithm 
as outlined in sec. 4.2. Given a system of linear equations, the algorithm of Solve 
can be described as follows: 

1. Calculate the ordering of the terms using FORM [5]. This involves calculating 
t, r and s for every integral according to their definition and counting the 
number of terms in each coefficient. Also, the number of terms in each equation 
is computed approximately. 

2. Find the integral with the highest t, r, s, with the smallest coefficient in the 
smallest equation. 

3. Solve its equation for this integral. The step 



involves inverting the coefficient cq. It is crucial to factor the coefficient if 
it contains symbols like d, the dimension of space-time. This is done using 
Maple [35]. 

4. Substitute the equation into the system. Store it for later use (see step 7). 

5. Perform partial fractioning to limit the growth of the size of the coefficients. 
Check that this removes the equation picked in step 2 from the system. This 
and the previous step are computationally expensive. 

6. Repeat all previous steps until no more equations are left. The system is then 
in a triagonal form and all identities are extracted. 

7. Substitute identities with some t, r and s that were extracted at a later cycle 
into those with the same t, r and s that were extracted earlier. 



4.6 The Code Solve 




(4.28) 
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8. Read out the identities, sorted according to the t, r and s of their left-hand 
sides; resubstitute them into the original system of IBPs to check that they 
indeed are a solution to it. 

Solve is written in C. It acts as a coordinator between various programs that 
are invoked to perform specific tasks; Maple for example is invoked to invert and 
factor the coefficient Cq. Then Solve writes FORM code for loading the system of 
equations, and Maple's output is used to write code to substitute the identity (4.28) 
into all others. For the actual substitution FORM is run on this code and Solve 
then processes FORM's output. Solve also keeps a file that stores the extracted 
identities and in the end rewrites them for the final output. 

The code Solve was used extensively in the preparation of this thesis. Solvers 
main virtue is the fact that it invokes FORM to do all algebraic manipulations. 
FORM is arguably the fastest computer algebra system to work with large symbolic 
expressions and so for large systems of IBP identities. Solve is quite efficient. The 
part of Solvers code that computes the ordering of the terms in the system of equa- 
tions is written as a plain FORM include file and thus can easily be adapted by the 
user. 



4.7 Improvements on Solve 

A number of extensions and improvements have been introduced while trying to 
make Solve fit the needs of the form factor reduction. Most of them try to address 
the fact that the system of linear equations for the three-loop form factor grows 
large in intermediate steps of the solution process, and thus Solve has to handle 
large expressions as efficiently as possible. These improvements are not specific to 
the form factor calculation and thus calling them improvements of Solve itself seems 
justified. In the following, we will comment on improvements of several steps in 
Solvers algorithm outlined above and we also introduce a few new features: 

l.,2. Computing t, r and s for a given topology T is actually quite fast. Still, there 
is no need to recompute them for all topologies in every step as suggested 
in Solvers algorithm. Rather, we may multiply cq by a new symbol and 
update t, r and s only for integrals that come with a factor A. For the next 
step, all A's are removed by setting them equal to one. 

Computing the number of terms of a given subexpression in FORM is quite 
challenging: For efficiency reasons, FORM has no general possibility to operate 
on several terms at the same time. Counting terms can be implemented using 
form's collect instruction to cast all terms to be counted into a single 
function's argument. Then it is possible to count the number of terms, e.g. 
using splitarg. However, FORM limits the size and number of arguments, 
and for long equations this method of counting terms does not seem viable. 
Earlier versions of Solve thus could not distinguish the size even of moderately 
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large equations. An alternative that turned out to be rather fast is based on 
the idea that we do not actually need the size of all equations; instead, we 
only need to identify among the equations with the smallest t, r, s and with 
the smallest coefficient the topology that is part of the smallest equation. 
Therefore we propose the following specific algorithm to perform the task in 
step 2 of Salve's algorithm: 

l.,2.-i) Compute t, r and s if needed. To count the number of terms in each 
coefficient, use bracket, collect and nterms to count the number of 
terms in the coefficient. If collect fails^, set the number of terms in the 
coefficient to a large value. 

l.,2.-n) To extract the smallest equation, take a copy of the system of IBP iden- 
tities and set all coefficients of all integrals equal to one. This drastically 
reduces the size of the system, making it fast to handle. 

l.,2.-iii) Apply a . sort. The coefficient of each integral is now equal to the num- 
ber of terms in the equation. Use bracket and collect to move this 
coefficient into the argument of the function according to which FORM 
sorts the terms. 

l.,2.-iv) Use bracket, and then f irstbracket_ to extract the label of the equa- 
tion with the smallest r, s, number of terms in the coefficient and 
number of terms in the equation. 

l.,2.-v) Use dollar variables^ to pass the information on the equation just selected 
to every term in the system of equations. Adapt the ordering information 
such that the given equation is sorted before all others. 

Thus, information that is present locally with every integral (i.e. t, r, s) 
is computed locally. Information that has to be collected from all over an 
equation is not computed for every equation, but instead only the term needed 
for the next substitution step is singled out. The information on the size of 
the coefficients is actually not local either and can be computed very much in 
the same way as the number of terms in the equation. But the difference to 
the method used above does not seem to be big because large coefficients tend 
not to be factorisable in a simple way. This provides a strong reason not to 
use them, since partial fractioning is very expensive in that case. 

3. When substituting in an identity, it is applied to the system containing the 
equation that gave rise to the identity. It is thus crucial that this equation 
indeed reads = after substitution and is thus removed from the system. 
Solve explicitly checks this after every substitution of an equation. This in 
turn requires and insures that partial fractioning is done properly after every 

^form's collect statement has a special second argument to deal with this case. 
^For earlier versions of FORM, one may also use multiply. 
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substitution cycle. 

Now our partial fractioning routine in FORM is not able to handle completely 
general denominators. Factoring I/cq into small factors that can be handled 
by Parfrac.prc might not be possible, and the above-mentioned check would 
fail. Thus a subroutine has been devised written using Mathematica's [34] 
powerful language; it has to 

a) try to factor the coefficient; 

b) analyse the factorised coefficient and decide what Solve should do with 
the current identity. 

The information on the outcome of the above operations is handed back to 
Solve, together with the inverted coefficient (in case the inversion was success- 
ful). Solve then writes code according to the following four scenarios: 

a) The coefficient is a pure number. No partial fractioning is needed in this 
case. 

b) The coefficient is indeed factored into terms that can and must be handled 
by Parfrac.prc. In this case do invoke the partial fractioning routine. 

c) The coefficient is not factored into terms that can be handled by Parf rac . pre. 
Write code that ehminates the equation from the system. For the final 
check of the solution of all IBPs, a list of dropped equations must be 
kept. 

d) The coefficient is vanishing. This may happen since partial fractioning is 
not used at every step, and FORM does not recognise terms like — 
1 — as vanishing. No substitution has to be made, but Parfrac.prc 
is invoked. 

4. Here, only minor changes were needed. The factor A was already introduced 
earlier. Also, FORM cannot substitute large identities using a single id state- 
ment. Thus, insted of 

id T=Tl+T2 + ... + T50+T51 + ...; 
. sort 

constructions like 

i d T=-fTl+T2 + . . . T5(HsumUp ; 
. sort 

id sumUp=sumUp-|-T51 -I- . . . ; 
. sort 

were needed already to cope with medium sized expressions. 
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5. For small systems of IBP identities, their solution finally contains only a rather 
limited set of factors in their coefficients. Their partial fractioning may be en- 
tered by hand. This approach revealed its weaknesses in the application of 
Solve to the three-loop QCD quark form factor. The growing number of dif- 
ferent factors appearing in the denominator of the solution of large systems 
made it difficult to anticipate all possible denominators. Rather, it became 
necessary to systematically perform partial fractioning for a general polyno- 
mial that contains factors of the form + cad, OL\^ai G Q. Instead of using 
symbols such as [o? — 4], [3 * o? -|- 2]~^ etc. and implement partial fractioning as 

s [d-4] ,[3*d + 2]; 

id [d-4]*[3*d + 2]"(-l) = l/3-14/3*[3*d + 2]"(-l); 

for all combinations, we switch to FORM's functions and use b(d,-4) and 
l/3*B(d,2/3) = 1/3 * 1/ (ci-F 2/3) instead. Then, partial fractioning may be 
implemented as 

cf b,B; 

id b(nl?)*B(nl? ,n2?)=l-b(n2)*B(nl ,n2) ; 

As can be seen from this example, it is crucial to enforce an ordering in the 
arguments of B(..): if nl is a number and n2 = rf, the above substitution should 
not be made- otherwise FORM runs into an endless loop. 

The use of FORM's repeat . . . endrepeat statement produced stack overflows 
when applied to large expressions and had to be replaced by 

#do k=l,l 

if ( match(<Pattern >) ); 

redefine k "0"; 

id <Pattern > = ..; 
endif ; 
. sort 
#enddo 

7. In the final check, identities that were dropped have to be accounted for. These 
equations are of course not solved by the identities collected, so the check will 
fail unless the dropped equations are removed separately. 

Taken together, the above improvements allowed for the solution of systems of 
IBPs that contained larger equations and more complicated denominators. Impor- 
tant improvements were the reliable and general partial fractioning of the coefficients 
of the integrals and the accurate and fast counting of the number of terms in large 
equations. 
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4.8 Extensions of Solve 

The code Solve as described in sections 4.6 and 4.7 is essentially an implementation 
of Gauss' algorithm to solve systems of linear equations. We have described several 
improvements on the efficiency, generality and stability. Still, further improvements 
are necessary to handle the size of the systems of IBP identities that were encoun- 
tered. Most importantly we have not yet made any use of our knowledge about the 
structure of IBP identities. The following facts are crucial: 

• Systems of IBP identities are sparse. 

• There are subsystems that have non-vanishing coefficients (almost) only in 
specific blocks in the coefficient matrix. 

The first point follows immediately from the number of IBP identities per topol- 
ogy and r-s-block (see tab. 4.2) and the fact that IBP identities are generated from 
seed integrals using differentiation, eq. (4.3). Using the chain and product rule, 
one can estimate that an IBP identity will typically contain O (10) terms. Tab. 4.2 
shows clearly that there are usually many more equations than that, and so most of 
the entries in the matrix of coefficients will be zero. 

The second point follows from the r-s structure of IBP identities: We have 
already classified integrals into t-r-s blocks, and we have seen in (4.11) that seed 
integrals from a given r-s-block mainly relate integrals of three r-s-blocks: Mostly 
r -|- 1, s -|- 1, but also r + l,s and r, s integrals. 

Considering these two facts, it seems natural to divide large systems of IBP 
identities into subsystems of identities that originate from seed integrals within a 
given r-,s-block. These identities mix strongly, and once their solutions are found, 
most of the substitutions of identities will be done. The overlap with the solution 
of a different r-s-block of seed integrals will be quite small; it is significant only for 
the r, s — 1, the r — 1, s and the r — 1, s — 1 block, see fig. 4.4. The next section thus 
presents routines that deal with the merging of such (sub-)systems of equations. 

We recall that dividing up a system of IBP identities into r-s-blocks is also useful 
in the light of sec. 4.4.1. There we saw that seed integrals in a higher r-s-block have 
more complicated 1-subtopologies and that their use should be avoided if possible. 
Thus it is in fact useful to merge and apply identities with an increasing number of 
r-s blocks of seed integrals, in contrast to simply merging all available r-s-blocks. 

4.8.1 Mergelds and Rewritelds 

The standard task that we would hke to describe in this section is the following: 
Given the solution of a system of IBP identities generated using seed integrals of a 
certain set of r-s blocks, and the solution of the IBP identities of an additional r-s 
block (usually adjacent in r-s-space); to find the solution of the combined system of 
IBP identities. 
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Figure 4.4: An illustration of how r-s-blocks of IBP identities interfere with each 
other. IBP identities generated from the r-s-block contain mainly integrals in the 
triangle (r, s), (r + 1, s), (r + 1, s + 1). They must be merged with identities stem- 
ming from three blocks with lower r and s: (r, s — 1), (r — 1, s) and (r — 1, s — 1). 
Identities from the r-s + 1-block also mix in (dashed line) but since their seed inte- 
grals have higher r and s, we may argue that IBP identities from r-s seed integrals 
mix into IBP identities from r-s -\- 1 seed integrals, and not vice versa. 



A simple way to achieve this is to read in the two systems of IBP identities and 
restart Solve. A routine that we call Mergelds scans the identities it receives and 
writes FORM code that builds up the FORM expression containing the system of 
equations: 

1 zeroes= sumlJp; 
. sort 

id sumUp=sumUp+Label ( 1 ) * ( Save (Ihsl) — Ihsl ); 
id lhsl=rhsl ; 
. sort 

id sumUp=sumUp+Label (2 ) * ( Save (lhs2) — lhs2); 
id Ihs2=rhs2; 

id Save(x?)=x; 
id sumUp=0; 

(the factor Label () distinguishes different equations). This system is then read into 
Solve that solves the system as described in section 4.6. 

There is one special situation where calling Solve is not necessary: Sometimes it 
is clear that it is necessary to substitute the identities of one set into the other, but 
not vice versa. This is the case for example when the left-hand sides of identities 
have different r or s. It is then enough to read in the identities of one set, apply 
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the second set of identities and write out the identities of the first set again. Solve 
is not needed in this case. We call this code Rewritelds: 

1 cf End; 

2 1 zeroes= 

3 + Label ( Ihsl ) 

4 + Label(lhs2) 

5 + . . . 

6 ; 

7 . sort 

8 id Label(lhs?) = Label(lhs)*(lhs+End); 

9 . sort 

10 ^include 1 stIdentitySet 

11 ^include 2 ndldentity Set 
i2#call Parfrac.prc 

13 .sort 

14 b Label ; 

15 pr int+s ; 

The contents of the bracket of Label (lhs)*( . . ) is the right-hand-side of the equa- 
tion, with the second set of identities substituted in. The last term of each right- 
hand-side is End, making it easy for a simple Perl script to read out the identities. 
We may note here that the generahsation to more that two sets of identities is 
straightforward. 

We would also like to point out that the two code samples above illustrate two 
different ways of reading in identities. The second one must be used with care, since 
if lhsl=lhs2, id lhsl=rhsl will be apphed both to Ihsl and lhs2 and the second 
identity is lost. Therefore, the above (outdated) version of Rewritelds only works on 
merged identities (where all left-hand sides are different) . Of course, this limitation 
is easily removed: 

1 cf End; 

2 1 zeroes= sumUp; 

3 . sort 

4 id sumUp=sumUp-|-Label ( Ihsl ) * ( lhsl-|-End ) ; 

5 id lhsl=rhsl ; 

6 . sort 

7 id sumUp=sumUp-|-Label ( lhs2 ) * ( Ihs2-|-End ) ; 

8 id Ihs2=rhs2 ; 

9 ... 

10 id sumUp=0; 

11 #call Parfrac.prc 

12 .sort 

13 b Label ; 

14 print-|-s ; 
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Mergelds and Rewritelds perform partial tasks of solving systems of equations: 
Mergelds substitutes into each other equations with the same left-hand-side using 

Solve. Rewritelds substitutes equations that reduce the integrals on the right-hand- 
side. Rewriting identities using Rewritelds often condenses the identities noticeably 
which is important for keeping small the size of terms involving subtopologies of 
subsequent (e.g. higher-t) IBP identities. 

4.8.2 PerlMerge 

Mergelds and Rewritelds are two routines that provide much flexibility in the choice 
of how to combine subsets of IBP identities. In our context, they are mainly useful for 
processing large systems of IBP identities that are split into r-s-blocks, which can be 
solved more easily. Using this flexibility however requires a lot of manual interaction, 
and in this section we will present a new routine, called PerlMerge that takes over 
some of this manual work. To illustrate the need for further automatisation, consider 
the reduction of the t — 7 topology with ID — 734: merging all identities from seed 
integrals with r = 7, s < 2 with those of r = 7, s = 3, we have the following t, r, 
s-tuples for the left-hand sides of the solutions: For r = 7, s < 2, 



and for r = 7, s = 3, 
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We clearly see that to merge these two sets of equations, there is no need to run 
Solve over the identities with r = 8, s = 4 again: There are not two sets of identities 
that would have to be merged. We may leave the r = 8, s = 4 identities aside 
and start merging the two r = 8, s = 3 sets of identities. There it happens that 
identities with r < 8 or s < 3 emerge. Those identities would have to be merged 
with other identities from the two sets. But in any case it is not necessary to merge 
all identities with r < 8, s < 3 at the same time. In this way Solve can be sped up 
significantly, since Solve has to process much smaller sets of identities at any given 
time: Identities whose Ihs have a smaller t-r-s than the ones currently processed do 
not need to be loaded into Solve yet. 
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The code (written in Perl) accepts a list of identity files as generated by Solve. 
Here is the algorithm: 

1. Extract the t, r and s of the left-hand sides of all identity files in the list. Sort 
the identity files accordingly 

2. Determine the highest t-r-s triple of all identity files in the list. 

3. If there are two or more files with this t-r-s triple, merge these identities. 
Otherwise copy the identity file to the directory storing the results. 

4. Update the file fist with the files containing the solutions of the merging process 
in the previous step. Update the list of t-r-s triples. 

5. Cycle through the fist of sorted t-r-s triples and determine the next lower 
triple. Go to step 3. 

The above algorithm is designed to be used with large systems of equations which 
is broken into smaller subsystems. They are then merged successively. As already 
mentioned one possible way of breaking up the systems is into r-s-blocks. 

One further application concerns the solution of a single (large) r-s-block e.g. 
t = 6, r = 8, s = 2 (or s = 3). It might help breaking it up into smaller, tractable 
parts again. Their solution can be parallelised. We then use PerlMerge to merge 
these subsystems gradually. When doing so, at some point the solution will start 
breaking up into parts with different r and s for the left-hand sides. They do not 
have to be treated at the same time anymore. Thus the systems to be merged with 
Perlmerge will be smaller than the original system. 

We found that for the partitionings we tried the gain was small. For example, 
dividing the system into the 16 pieces that arise when each IBP or LI operator 
in (4.3) is applied to all seed integrals, the systems typically decrease in size by 
about 10%. At some point the merged subsystems of equations again become large, 
and the gain of ~ 10% did not help much there. 

We should stress again that in the cases where it was possible to obtain the 
solutions of the single r-s blocks of identities, PerlMerge and Rewritelds did speed 
up the solutions of those topologies significantly, to the point that e.g. the solution 
for t = 7, r<8,s<2 (without having the necessary reduction of all t = 6 integrals 
at hand) would have taken much longer (if it could be obtained at all) using just a 
single run of Mergelds. 

Fig. 4.5 shows a fiow chart summarising how we approached the solution of large 
systems of IBP identities, combining all the tools discussed in this thesis. 
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Figure 4.5: A flow-chart of liow PerlMerge, Solve and Rewritelds work to- 
getlier on large systems of IBP identities. The system is split up into r-s-blocks 
(ri, si) , . . . , (r„, Sn)- A given r-s block's IBP identities are solved and then merged 
with lower-r-s blocks. This yields identities whose left-hand-sides (Ihs) are all differ- 
ent. Rewritelds then recursively substitutes these identities into each others right- 
hand-side (rhs) to yield the most compact expression possible. Then the next r-s- 
block is solved and added to the system. 
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Chapter 5 



Reduction 
Topologies 



of the 5-Propagator 



In the previous chapters we have set up a detailed analysis of how the Laporta 
algorithm could be employed for the reduction of the three-loop quark form factor 
to master integrals. We saw that this task is quite involved due to the large number 
of subtopologies and IBP identities that are present at the level of topologies with 
7 or less propagators. Quite a variety of tools were developed in order to tackle 
different types of topologies: Reducible topologies may be reduced using symbolic 
IBP identities, some topologies turned out to be integrable loop-by-loop, and for 
irreducible topologies that are not integrable loop-by-loop we devised routines such 
as PerlMerge, Mergelds and Rewritelds that all aim to tunc up Solve in the hope to 
be able to solve the large systems of IBP identities needed to reduce the form factor 
to master integrals. 

Here we should stress that for the topologies (at t = 6 and t = 7) where explicit 
numeric IBP identities have to be solved and where the number of IBP identities 
grows large even for singcl r-s-blocks of seed integrals, it is crucial that (lower-t) 
subtopologies are reduced as far as possible before starting to solve the system of 
IBP identities. During the solution of a system of IBP identities, many identities are 
substituted into others. If all these identities contain many lower t subtopologies, the 
equations quickly become very large. It is thus mandatory to reduce subtopologies 
to their master integrals first in order to simplify the part of the identity containing 
lower-t subtopologies. 

As we noted in sec. 4.3, at the level of t = 5 all topologies can be integrated loop- 
by-loop. At the same time the number of seed integrals is so large that using IBP 
identities is clearly hopeless. This chapter is thus going to describe how t — 5 {or 
t — A) topologies can be calculated loop-by-loop. We first discuss how to evaluate a 
single such integral in sec. 5.1. Since we need a lot of these integrals, we then also 
discuss the automatisation of their evaluation in sec. 5.3. To keep the size of the 
results small, we will rewrite all integrals in terms of their r = t, s = integral. 
These might not be master integrals. In section 5.4 we will give their reduction 
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identities to master integrals. 



5.1 Integration of Multi-Loop Bubble-Insertion Ver- 
tex Integrals 

We start with the following basic results for the 1-loop bubble and triangle integral: 

'"^^^^ = IWf{-k^-^r|r{-ik-if-^rJf ^'''^ 
= {-f-iv)'-''-'fa,, (5.2) 

^"^"^^^ ^ / (2n)' (-fc2 _ ,^)" (_ ^k-p,f -ir^f{-{k- qf -iT^f ^^'^^ 
- {-<l'-^v)'^-''-'-'9aA. (5-4) 

where a small imaginary part i 77 is understood in every scalar product or propagator. 
fa,i3 and 5'a,/3,7 consist of Euler gamma functions: 

f ^ z r(f-a)r(f-/3)r(-f + a + /3) 

^"-'^ (4vr)'^/2 V{a)T{d-a-(3)T{(3) ^^"^^ 

^ ^ r(f-a-/3)r(f-/3-7)r(-f + a + /3 + 7) 

(47r)''/^ r (a) r (d - a - /3 - 7) r (7) ■ ^'■''^ 

Eq. (5.2) replaces the two propagators of a bubble insertion, and \/ [k — i)"^ 
with a single propagator that has the same external momentum flow k—{k — i) = i. 
Thus for integrals that consist of bubble insertions into the 1-loop two- or three-point 
function (such as fig. 5.1), the number of loops may be reduced by replacing repeat- 
edly the left-hand-side of (5.2) by its right-hand side. The propagators that replace 
the bubble insertion then have an exponent that depends on the dimensionality d, 
but do not otherwise complicate the integral. 

From this follows a straightforward way to implement the evaluation of multi- 
loop integrals that consist entirely of two-point bubble insertions into a two- or 
three-point function: 

1. Count the number of times ki, k2 and k-i, appear in the integrand. Choose as the 
next loop momentum one that appears in a minimal number of propagators^. 
Call it k'. 



^The way our loop momenta are assigned to the propagators is crucial for this way of singling 
out the first two-point function to be integrated. Without tliis, one would have to consider the 
topological information of the diagram to find two-point insertions (i.e. two propagators connecting 
the same two vertices). 
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Figure 5.1: A sample 5-propagator topology that can be integrated loop-by-loop. 
After integrating out two bubble-insertions (loop momenta k2 and ks) , the standard 
vertex integral (with non-integer exponents) is recovered. 



2. Shift the loop momentum k' ^ k to fit the propagators into the standard form 
(5.2) or (5.4). Determine £. 

3. Collect the scalar products in the numerator that depend on k, and isolate the 
tensor loop integral (J d'^-kk • Pi — > Pi / cfikk^, etc.). 

4. Perform the tensor reduction, as will be described below in sec. 5.2. This 
yields the tensor integral as a scalar integral times a tensor formed out of I in 
eq. (5.2) or pi and q in (5.4). 

5. Evaluate the loop integral. This yields a new propagator (— £^ —irf)'^ " ^ 
times a coefficient fa^p or ga,^,^- 

6. Contract the tensor in i with the rest of the amphtude that has been split 
apart in step 3. Try to cancel scalar products against propagators. 

7. Go to step 1, proceed with the evaluation of the next loop. 

To illustrate the above steps (without tensor reduction) , consider as an example 
the diagram shown in fig. 5.1. We proceed as follows: ki appears in three different 
propagators, but k2 and k^ only in two and thus we clearly have two two-point 
functions to integrate. As a first step, either the integration over k2 or k^ is possible. 
We choose to integrate over k2- The propagators containing k2, 1/ (/c2 — Pi)^ and 
1/ (^2 — ks)'^ have to be shifted e.g. by k2 — > /c2 +Pi to fall into the form (5.2), with 
k = k2 and i = k-^—pi . Note that is indeed a member of the propagator set (3.9). 
Thus in the next iteration we are able to integrate over k^. The exponent of one of 
the propagators is now d-dependent. Finally, in a third iteration we encounter the 
integral (5.4). 
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For fig. 5.1 this procedure yields two fa,i3 coefficients and one ga,i3,'y Of course 
the dependence on in our massless calculation is fixed by the dimension of the 
integrals: There will always be a factor (— — ir])^^'^ 

The functions fa,i3 and 5'a,/3,7 in (5.5) and (5.6) are rather lengthy. But since only 
Gamma functions are involved, these expressions become quite short when divided 
by /i,i and gi,iX- Using the Pochhammer symbol 

^"^"^^T(^^^^" + ^^---^" + ''"^^ ' ^^-^^ 

we have for example 

fa,l3 _ 1 (2-f),+^_2((i-«-/3)„+^_2 



A,, r(a)r(/3) 



(5.8) 



Equivalently, we may divide any integral that can be evaluated loop-by-loop by 
its scalar r = t integral. Using partial fractioning then reduces the expression to a 
form similar to those obtained by IBP identities: rational polynomials in d times 
the scalar r — t integral. 



5.2 Tensor Reduction 

For integrals with s < 4 the tensors that might appear in the above procedure are 

Uj hj ^ J.^ rv ^ rv rv ^ rv rv rv rv ro nj r 



/ 



{2nf (_A;2)"(-(A;-£)2)^ 
d'^k {l,k^',k^'k'',k^'k''kP,k^'k''kPk''} 



(5.9) 



(5.10) 



In this section we derive the reduction of these integrals to the scalar integrals (5.2) 
and (5.4). As an explicit example, we take the integral 

d'^k k^k" 



but allow for a natural generalisation to higher tensors. We first construct a list 
containing all possible Lorentz 2-tensors that might arise in this integral: 

V^'' = ( ) • (5-12) 
The tensor integral can now be written as 

i!::p^v^'^rc, (5.13) 
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where C = ^ ^ is the coefficient vector that is to be found. Multiplying (5.13) 



from the left with V^^, and inverting, we get 

Substituting this back into (5.13) we arrive at 

= V^'^''. {V^-.V^^) ' .Ka./S . (5.15) 

The result (5.15) holds not just for general a and f3, but generalises without difficulty 
to higher rank tensors. However, the matrix {V'^^'.V^^ has entries that depend on 
the dimensionality d = g^ud'^", and thus its inversion becomes cumbersome for 
tensors of rank higher than four. 
As an example, we get 

^a,/3 (?) = Ha-l,p (?) + Ia,P-l (?) + ?'/a,/3 (?)) ■ (5.16) 

The result for I^'^^ is 

+ Ia-2,P (?) - 2/a-l,/3-l (?) + 4,^-2 (?)^ ■ 

5.3 Parallelising the Evaluation of 5-Propagator 
Integrals 

So far we have discussed the evaluation of a single t = 5 integral. We will next 
comment on how to organise the calculation of all t = 5 integrals needed for the 
reduction and how to store the results. 

The 5-propagator integrals appear as subtopologies in the t — 6 IBP identities. 
If these would be solved without a prior reduction of the 5-propagator topologies, 
the solutions would contain unmanageably large expressions involving various 5- 
propagator integrals. Thus in order to keep those systems small, it is mandatory 
to reduce the t — 5 integrals first. Some of the t — 5 topologies are reducible, in 
which case they reduce (using symbolic IBP identities) to the only non-vanishing 
t = 4 (irreducible) topology. The irreducible t = 5 and t = 4 topologies are shown 
in fig. 5.2. Their integrals must be evaluated as described above. 

For the evaluation of integrable subtopologies there are two possible approaches: 
Either we extract the integrals needed from the t — 6 system of IBP identities and 
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Figure 5.2: The 5- and 4-propagator irreducible subtopologies. All t < 5 integrals 
can be integrated loop-by-loop and ultimately be expressed as a linear combination 
of these five master integrals. This is the form most suited for a substitution into 
the t — 6 system of IBP identities. 



integrate those topologies only once the t = 6 system is set up. This results in a 
possibly rather small set of integrals that has to be computed for a given set of IBP 
identities; but there is a danger of recomputing the same integral several times if 
they show up in different sets of IBP identities. A second possibility is to tabulate 
all integrals in a given r-s range (say i<r<9, 0<s<4) and compute them all, 
but only once. 

The choice depends on the context. For loop-by- loop integrable t = 7,8 topolo- 
gies^, it was affordable not to integrate them when they appear as subtopologies in 
IBP identities of higher-i topologies. There we would first reduce the integrals and 
integrate the (less than 1000) loop-by-loop integrable integrals only once they show 
up in the reduction result for the form factor. An important side remark here is that 
loop-by-loop integrable topologies may not only be irreducible themselves, but they 
may also be reducible and thus contain difficult (irreducible and not loop-by-loop 
integrable) subtopologies. Therefore it is important to detect these topologies in 
the result of the form factor and integrate them instead of reducing them further. 
This of course requires determining for every topology whether it is integrable or 
not (which we did by hand). 

For the 5 and 6 propagator loop-by-loop integrable topologies we chose instead to 
evaluate all r < 9, s < 4 integrals of the irreducible t = 4, 5 topologies needed. This 
job can be nicely parallelised and the computational effort is acceptable: Table 4.3 
counts a total number of 62370 integrals per topology with r < 9, s < 4. Using 



^There is no loop-by-loop integrable t = 9 topology needed for the quark form factor. 



5.4. REDUCTION EQUATIONS FOR T = 5, 6 INTEGRALS 
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table 4.1, we are to compute 6 x 62370 = 374220 integrals. The computation takes 
at most a few seconds per integral, so we need at most a few CPU-weeks. The results 
are then compiled into several FORM tablebases [36]. This has the advantage of a 
short look-up time and the results can easily be reused whenever needed. Also, the 
tablebases contain the results in a compressed form. The total size of the tablebases 
was 221 MB. 

For the implementation of the algorithm of sec. 5.1, one important restriction 
to keep in mind is that the computation should not be implemented in Mathemat- 
ica [34] or any other software with expensive licenses, since they do not allow for 
parallelisation. Here again FORM [5] was our product of choice. 

5.4 Reduction Equations for Scalar t = 5,6 Inte- 
grals to Master Integrals 

In this chapter we have so far discussed the evaluation of (mainly t = 5 or t = 6) 
integrals that are integrable loop- by-loop. After their evaluation in terms of Euler 
gamma functions, we used eq. (5.8) to express them in terms of their r — t, s — 
integral. In this way all gamma functions cancelled, and the coefficient consisted of 
a polynomial in d, eq. (5.8), similarly to IBP identities. The r = t, s = integral 
however might still be reducible. This can easily be found out by solving a few IBP 
identities. We will next give the reduction equations for the loop-by-loop integrable 
topologies that appeared as subtopologies oi t — 6 and t — 7 topologies. The t — 5 
reducible subtopologies obviously reduce to 7^72 = B4: 



For the t = 6 reducible subtopologies, the following master integrals appear on the 
right-hand-side: Tigg = B^a, 1^2 = ^5,2, ^33 = ^5,1, ^662 = C'2,3- 
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Chapter 6 

Current Status and Outlook 



In chapters 2 to 5 we have set up tools that we hoped to be sufficient to reduce 
the three-loop quark form factor to master integrals. We will next describe our 
experience made so far with the reduction of the three-loop quark form factor. Since 
we have not yet completed the reduction, sec. 6.1 will outline the progress made so 
far and highlight in which form the difficulties anticipated above were encountered. 
Sec. 6.2 will give an outlook on possible future applications of our tools. In sec. 6.3 
we present a partial result. 

6.1 Applying Solve to the Three-Loop Form Fac- 
tor 

As already stated in sec. 4.4, the reduction of t = 9 and t = 8 topologies do not pose 
a problem, though the improvements in Solve did help (for example, large, non- 
factorisable coefficients are already encountered in the reduction oit — 8 topologies 
and would have caused earlier versions of Solve to interrupt). It was affordable to 
take a top-down approach, where no t—1 subtopologies are reduced before solving the 
system from t-propagator seed integrals. Not reducing the 7-propagator subtopolo- 
gies has the disadvantage that the intermediate expressions in the reduction of t = 8 
topologies become rather lengthy. But in this way the following information could 
be gained: At the level oi t — 7 propagator topologies, no r = 9, s = 4 topologies 
were left. In fact, a single r = 9, s = 3 topology remained that was already present 
in the amplitude for the form factor. It turned out to be reducible using symbolic 
IBP identities, and not = 6, r = 9,s = 3 subtopologies emerged. In fact, after 
applying also the symbolic reduction equations (and integrating loop-by-loop where 
possible) no t — 7, r — 9 integrals were left. The worst irreducible topologies that 
could not be integrated loop-by-loop had (t — 7) r = 8, s = 3. 

One might have hoped that at t = 7, using seed integrals up to r = 7, s = 3 
only, all such t = 7,r = 8, s = 3 integrals would reduce. This was not the case, so 
seed integrals with t — 7,r — 8, s — 3 must be solved. 
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At this point we tried to start the proper bottom- up approach: 1-subtopologies 
should be reduced to a form as compact as possible and substituted into the t- 
propagator topologies. Wc have already seen that 5-propagator topologies can sim- 
ply be integrated out and then related to master integrals. Ideally the next step 
would be to reduce all needed t = 6 topologies (11 in total) up to r = 8, s = 3, such 
that the t — 7 IBP identities contained the least possible number oit <7 integrals. 
This so far was not possible because of the large number oi t — 6, r — S, s — 3 
IBP identities that have to be solved (as was illustrated in tab. 4.2). Indeed, even 
identities from t = 6,r = 7,s = A seed integrals sometimes grow too large for an 
efficient handling using Solve as described. 

This is where the reduction stands at the moment. The amplitude still contains 
O (10) t — 7,r — S, s — 3 integrals, and presumably once they are reduced a few 
complicated t = 6, r = 8, s = 3 integrals will be present. Also t = 6, r = 8, 
s = integrals that are not reduced using t — 6, r — 8, s — seed integrals might 
eventually pose a problem. 

An idea still to be investigated is that equations with a long tail of t — 1 subtopolo- 
gies could be shortened by introducing abbreviations. The tools to do so are already 
set up, and first experiences are positive: With some patience, the reduction of IBP 
identities from t = 7,r = 8, s = 3 seed integrals might become possible without 
reducing any t = 6 subtopologies by simply hiding them in (nested) abbreviations. 
Also, these systems seem to contain enough equations to overconstrain the whole 
t — 7,r — 8, s — 3 block of integrals, such that once a solution could be obtained, 
it would reduce all these integrals instead of just relating them to each other. Using 
the solutions of lower-r-s seed integrals, also the remaining (e.g. t — 7,r — 8, s — 2) 
integrals might find a reduction equation. 

6.2 Outlook 

So far we have dealt exclusively with one out of three auxiliary topologies needed 
in the reduction of the three-loop quark form factor. Once the reduction of the 
planar auxiliary topology to master integrals is achieved the reduction of the re- 
maining two auxiliary topologies is clearly the most immediate next step necessary 
to obtain a complete (gauge invariant) result of the three-loop quark form factor. 
We have pointed out however that this task will be a smaller one since most of the 
subtopologies encountered there will be the planar ones discussed here already. 

Upon completion of this reduction, several immediate extensions are possible 
using the methods devised in this thesis. They can be applied without significant 
changes to any massless three-loop QCD vertex correction. The three-loop QCD 
form factors of the gluon, Higgs- and Z-boson are thus of primary interest. These 
calculations will profit heavily from the results of the photon form factor. In the 
case of the gluon form factor new topologies are expected to arise. The Higgs form 
factor will probably contain an additional scalar product in the numerator of its 
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t — 9 integrals which might pose a formidable challenge for the reduction. On the 
other hand the Z-boson form factor is expected to be very similar to the photon 
form factor since the difference only lies in the vertex Feynman rules and does not 
affect the topologies of the diagrams. 

More generally speaking, in the quest for higher order perturbative calculations 
the niche of the tools discussed in this thesis lies in processes where the systems 
of IBP identities are too large for other implementations of the Laporta algorithm 
such as AIR [37] . AIR is very user- friendly and well tested but the solution of the 
system of IBP identities is carried out entirely in MAPLE and not in FORM. For 
large system this is expected to be a problem. 

Of course, Solve and its extensions devised in this thesis do not rely on a par- 
ticular physical process. If the Qgraf Feynman rules are implemented, non-QCD 
(e.g. electroweak, supersymmetric,...) processes can be considered, too. Occasion- 
ally we have specialised in vertex graphs, e.g. in the the standardisation of the loop 
integrals. But the reduction itself only involves the solution of a general system of 
linear equations. Setting up and solving the system of equations only involves minor 
adjustments (e.g. of the coefficients of the momenta that define the propagators of 
the auxiliary topology). The reduction is then automated to a large degree but a 
lot of control over the details of the reduction can be maintained if necessary. 

Summarizing we might say that the Laporta algorithm has brought about great 
advantages for many applications, but it also shows its limits: The prize for the 
automatization of the reduction process is the need to solve a large system of linear 
equations. These systems quickly become large for present-day computers. Still, we 
are not aware of any other method with which an analytic result for three-loop form 
factors would be in reach. 



6.3 Partial Results 

The reduction of the 3-loop quark form factor is yet to be completed. However, we 
may already present one partial result. Collaborators reported that the evaluation of 
the 9-propagator master integral (with topology ID=1790) is actually much harder 
than evaluating the integral 

Ti79o (0, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 0) , (6.1) 

i.e. the integral where the integrand has an additional scalar product (/c2 — Pi)^ in 
the numerator. The reduction of this integral is much easier than the reduction of 
the full amplitude for the form factor: At t = 6 the most difficult integrals that need 
to be reduced have r = 8 and s = 0. The corresponding reduction equations could 
be found, and we successfully cross-checked our result with our collaborators who 
obtained the same result independently using AIR [37]. We obtain 
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6.4 Summary and Conclusions 

In this thesis we have attempted to calculate a part of the three-loop QCD contribu- 
tion to the process 7* qq. We generated the 244 Feynman diagrams and projected 
them onto the form factor. The integrals in this expression had to be analysed and 
standardised in a semi- automated way, since even after standardisation the expres- 
sion contained 2951 complicated three-loop integrals. The ultimate goal was to use 
integration-by-parts identities -linear equations in the unknown integrals- that in 
principle relate all integrals to a set of master integrals. Most of these master inte- 
grals have already been computed by our collaborators [29]. But the reduction to 
master integrals proved difficult due to the large number of equations to be solved 
for the irreducible topologies. A final result is still missing, and so the main chap- 
ters in this thesis described the various specialised methods and algorithms that we 
hoped to push to a sophistication level sufficiently high in order to achieve a result. 

In the meantime -after having finished this thesis- the quark and gluon form 
factors have been computed to three loops [38] . Congratulations to the authors! 
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Appendix A 
Qgraf 



Qgraf [4] is a tool for the fast generation of Feynman diagrams. After reading in 
the Feynman rules, it essentially provides the information that defines each diagram. 
This information is then printed according to a user-defined style file. This section 
describes in some (but not every) detail the files and configurations used for the 
generation of the three-loop quark form factor. 

Firstly, Qgraf needs a model file listing propagators and the possible vertices of 
the particles in the Lagrangian: 

Listing A.l: the file qgraf.dat 

I % propagators 

2 

3 [Q,Qbar,— ] 

4 [ Ghost , Antighost , —] 

5 [ Gluon , Gluon , -|-] 

6 [ Photon , Photon ,-|- , external] 

7 

8 % vertices 

9 

10 [ Qbar , Photon ,Q] 

II [ Qbar , Gluon ,Q] 

12 [ Gluon , Gluon , Gluon ] 

13 [ Gluon , Gluon , Gluon , Gluon ] 

14 [ Antighost , Gluon , Ghost ] 

This is a model file for QCD, plus an additional photon that couples to quark- 
antiquark pairs. The kej^word external in line 6 specifies that photons appear only 
as initial or final states in the diagram. Qgraf generates the diagrams using the 
Feynman rules in this model file. 

Qgraf output is user-defined according to the style file. The one we used reads 

Listing A. 2: the style file mystyle.sty for Qgraf 
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1 <prologue> 

2 * 

3 * file generated by <program> 

4 * 

5 <command_loopXcommand_line_loop >* <command_dataxback> 

6 <endxend>* 

7 

8 <diagram> 

9 

10 id sumup = sumup + 

11 Label(<diagram_index >)* 

12 (<signXsymmetry_factor >)* 

13 <in_loop> 

14 pol<field >(< field -index >,<momentum>)*<end> 

15 <out_loop> 

16 pol<field >(< field -index >,<momentum>)*<end> 

17 

18 <propagator_loop > 

19 prop<ficld >(<field_index >,<dual — field _index >, 

20 <backXmomentum>)* 

21 <backxend> 

22 <vertex_loop > 

23 vert <vertex_degree >(<vertex_index >, 

24 <backXray_loopXf ield >(<field_index >, 

25 <back><momentum>),<endxback> 

26 <back>)* 

27 <backxendxback >; 

28 .sort 

29 <epilogue > 

30 

31 <exit > 



In the diagram section (lines 8-28), Qgraf loops over all diagrams, external par- 
ticles, propagators and vertices and provides information for various quantities that 
define parts of the diagram. For example, in the propagator loop (lines 18-21), the 
quantities <f ield_index> and <dual-f ield_index> are provided for every propa- 
gator. They indicate which two fields are connected by a propagator. In the vertex 
and ray loop (lines 22-27), quantities like vertex_index and again <f ield_index> 
are available to define the vertices. 

Before being able to run Qgraf, the process must be specified for which the 
diagrams should be generated. This is done in the file qgraf.dat. The information 
given there includes the name of the output file (e.g. gammaqqbarSLoop . out), the 
style and model file to be used, the initial and final states of the process together 
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Figure A.l: The first diagram in the Qgraf output. The numbers associated to the 
propagators or vertices correspond to the (dual-)field numbers that Qgraf provides. 
They are used in the actual Qgraf output that is given in the hsting A.4. Using 
such information, the Feynman amplitude can be constructed. 

with their momenta, the numbers of loops, the root name to be used for the loop 
momentam and various options. 

Listing A. 3: the file qgraf.dat 

output = ' gammaqqbar3Loop . out ' ; 

style = ' mystyle . sty ' ; 

model = 'QCD_QED' ; 

in =Q[pl], Qbar[p2], Photon [mq]; 

out = ; 

loops = 3 ; 

loop_momentum = k ; 

options = onshell , nosnail ; 

Running Qgraf then produces the information that there are 244 diagrams. They 
are written down in the output file, a plain text file: 

Listing A.4: excerpt of the Qgraf output file gammaqqbar3Loop.out 



id sumup = sumup + 
Label ( 1) * 
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4 (+1/6)* 

6 

6 polQ(-l,pl)* 

7 polQbar ( — 3 ,p2 ) * 

8 polPhoton ( — 5 ,mq) * 

9 
10 
11 

12 propQ (2 , 1 , — kl ) * 

13 propQ (4 ,3 , — kl+mq) * 

14 propGluon (6 ,5 , — kl— pi ) * 

15 propGluon (8 ,7 , kl+pl ) * 

16 propGluon (10, 9, -kl+k2+k3-pl)* 

17 propGluon(12,ll,-k2)* 

18 propGluon (14 , 13 ,— k3) * 

19 

20 vert3 (1 ,Qbar(4 ,kl— mq) , Photon( — 5 ,mq) ,Q(1, — kl) )* 

21 vert 3 (2 ,Qbar(2 ,kl) , Gluon(5,-kl-pl ) ,Q(-l,pl) )* 

22 vert 3 (3 , Qbar ( — 3 ,p2 ) , Gluon (7 , kl+pl ) ,Q(3, — kl+mq) ) * 

23 vert4 (4, Gluon (6 , kl+pl) , Gluon (9 , -kl+k2+k3-pl ) , Gluon (11 , -k2 ) , 

Gluon(13,-k3) )* 

24 vert4 (5 ,Gluon(8,-kl-pl) , Gluon ( 10 , kl-k2-k3+pl ) , Gluon ( 12 , k2 ) , 

Gluon(14,k3)) ; 

26 . sort 

26 

27 id sumup = sumup + 

28 Label (2)* 

29 ... 

This diagram (the first in Qgraf^s output) is shown in fig. A.l. The momentum 
flow on vertices is always inwards, and the propagator propQ (4,3, -kl+mq) has the 
momentum —ki — q flowing from the field index 4 to the dual field index 3. Also 
note the symmetry factor and fcrmion loop sign in line 4. 

The above code is adapted to be read into a FORM file simply using 

1 1 expr=sumup ; 

2 ^include gammaqqbar3Loop . out 
The next steps are then: 

• Provide commuting factors such as 1/p^ for quark or gluon propagators, colour 
and electric charges for fermion vertices, 3- and 4-gluon-vertices. 

• Reconstruct external fermion lines: To this end find the vertex that con- 
tains the same field index (—3) as the antiquark (in our example: vertex 3 



in line 22). Then find the quark propagator that connects to the same vertex 
(propQ (4,3), -kl+mq) , fine 13) , then the vertex that connects to the next field 
index 4 (vertex 1, fine 20), etc., until the field index of the external quark (-1) 
is reached. 

• Reconstruct closed fermion loops similarly. 

• Project onto the form factor; evaluate Dirac traces and colour factors. 

• Collect the integrals into T's as defined in eq. (3.7). 

Flavour indices on propagators and vertices must be taken into account properly. If 

after projection onto the form factor in (2.3) the two Dirac matrices in the photon 
vertices 7^^, 7'' lie in the same Dirac trace, the charge factor is (e^), otherwise it 
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As described in sec. 3.4, quite a number of topologies are usually equivalent, since 
permutations of (loop or external) momenta change an integral's ID number, but not 
its value. The ID number defines the topology^ according to the auxihary topology 1 
in tab. 3.1. It is restricted to be smaller than 2^^ = 4096. Here we list the set of 
identifications that we obtained. Topologies that are not listed are either vanishing 
(e.g. all t < 3 integrals) or have t > 9 and thus are not present in the three-loop 
quark form factor. 

t = 4 

178,540,561,2074,2089 ^ 172 





677,2190,2211,2573,2579 ^ 662 

626,1193,1562,2156,2332 ^ 433 

604,1202,1585,2138,2345 ^ 428 

436,620,628,1178,1194,1561,1577,2154,2162,2329,2353 ^ 412 
242,434,625,796,817,1196,1564,2153,2330,3098,3113 ^ 236 
188,569,572,2105,2106 ^ 186 

182,541,565,668,684,2075,2091,2202,2226.2601,2609 ^ 174 
179,542,563,689,690,2078,2093,2217,2220,2586,2588 ^ 173 
557,692,2099,2218,2585 ^ 158 



^We have seen that several graphs may have the same ID number. The graph given in fig. A.l 
e.g. has two propagators with the same momentum. This is not reflected in its ID number. More 
precisely, the ID number only distinguishes r = t, s = integrals. 
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1699,2445,2469,2638,2646 ^ 1683 

1457,1626,1650,2396,2412 ^ 1449 

1642,2420 ^ 1433 

918,2275,2835,3214,3597 ^ 741 

933,1686,1701,2254,2446,2467,2637,2643,2829,3235,3603 ^ 726 

2221,2590 691 

685,693,694,2206,2219,2222,2227,2587,2589,2605,2611 ^ 670 

686,2203,2230,2603,2613 ^ 669 

679,2191,2215,2575,2583,2701,2702,2707,2710,2723,2725 ^ 663 

884,1258,1434,1436,1641,1644,1817,2417,2418,3178,3353 ^ 500 

881,1260,1820,3177,3354 ^ 498 

882,1257,1818,3180,3356 ^ 497 

860,1266,1452,1458,1628,1649,1841,2394,2409,3162,3369 ^ 492 

876,1242,1450,1460,1625,1652,1833,2393,2410,3186,3377 ^ 476 

636,1210,1593,2170,2361 ^ 444 

634,1209,1594,2172,2364 ^ 441 

629,748,924,1198,1565,2155,2290,2331,2865,3226,3625 ^ 438 

630,1195,1563,1690,1705,2158,2333,2460,2481,2668,2674 437 

627,754,945,1197,1566,2157,2284,2334,2844,3241,3610 ^ 435 

605,732,940,1206,1589,2139,2266,2347,2857,3250,3633 ^ 430 

606,1203,1587,1713,1714,2142,2349,2473,2476,2650,2652 ^ 429 

621,756,948,1182,1581,2163,2282,2355,2841,3242,3609 ^ 414 

622,1179,1579,1689,1706,2166,2357,2457,2484,2666,2676 ^ 413 

252,442,633,825,828,1212,1596,2169,2362,3129,3130 ^ 250 

246,797,821,1692,1708,2458,2482,2665,2673,3099,3115 ^ 238 

243,753,798,819,946,2281,2842,3102,3117,3244,3612 ^ 237 

813,1716,2474,2649,3123 222 

573,700,2107,2234,2617 ^ 190 

189,571,574,697,698,2109,2110,2233,2236,2618,2620 ^ 187 

183,543,567,2079,2095,2714,2716,2729,2732,2737,2738 ^ 175 

559,2103,2713,2730,2740 ^ 159 
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t^7 

2951,3719 ^ 2759 

2711,2727 2703 
1939,2533,2902,3469,3662 ^ 1763 

1957,2510,2893,3491,3667 ^ 1750 

1955,2509,2894,3493,3670 ^ 1747 

2477,2654 ^ 1715 

1707,2461,2485,2670,2678 ^ 1691 

1703,2447,2471,2639,2647,2766,2774,2957,2981,3731,3747 ^ 1687 
1522,1884,1905,3418,3433 ^ 1516 

1521,1882,1906,3420,3436 ^ 1513 

1500,1514,1524,1881,1897,1900,1908,3417,3434,3441,3442 ^ 1498 
1898,3444 ^ 1497 

1658,2428 ^ 1465 

1459,1630,1651,1778,1969,2398,2413,2540,2908,3497,3674 ^ 1453 
1461,1627,1654,1754,1961,2397,2414,2524,2924,3505,3698 ^ 1451 
1437,1643,1646,1770,1945,2421,2422,2548,2932,3481,3690 ^ 1435 
1438,1645,2419,3306,3865 ^ 1012 

1004,1454,1462,1629,1653,2395,2411,3290,3314,3881,3889 ^ 988 
997,1765,1942,2531,2899,3278,3299,3470,3661,3853,3859 ^ 982 
949,1694,1709,2286,2462,2483,2669,2675,2845,3243,3611 ^ 758 
947,2285,2846,3245,3614 ^ 755 

925,1693,1710,2294,2459,2486,2667,2677,2869,3227,3627 750 
757,926,950,2283,2291,2843,2867,3230,3246,3613,3629 749 
919,2279,2787,2789,2839,2963,2966,3215,3599,3725,3726 ^ 743 
941,1717,1718,2270,2475,2478,2651,2653,2861,3251,3635 ^ 734 
942,2267,2859,3254,3637 ^ 733 

935,2255,2765,2771,2831,2958,2979,3239,3607,3734,3749 ^ 727 
702,2235,2238,2619,2621 ^ 701 

2237,2622 ^ 699 

2223,2591,2718,2733,2739 ^ 695 

687,2207,2231,2607,2615,2715,2717,2731,2734,2741,2742 ^ 671 
892,1274,1466,1468,1657,1660,1849,2425,2426,3194,3385 ^ 508 
889,1276,1852,3193,3386 ^ 506 

890,1273,1850,3196,3388 ^ 505 

885,1262,1772,1821,1948,2546,2929,3179,3355,3482,3689 ^ 502 
886,1259,1769,1819,1946,2545,2930,3182,3357,3484,3692 ^ 501 
883,1009,1010,1261,1822,3181,3305,3308,3358,3866,3868 ^ 499 
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t — 7 (continued) 
861,1270,1756,1845,1964,2522,2921,3163,3371,3506,3697 ^ 494 

862,1267,1777,1843,1970,2537,2906,3166,3373,3500,3676 ^ 493 
877,1246,1780,1837,1972,2538,2905,3187,3379,3498,3673 ^ 478 
878,1243,1753,1835,1962,2521,2922,3190,3381,3508,3700 ^ 477 

637.764.956.1214.1597.2171.2298.2363.2873.3258.3641 ^ 446 
638,1211,1595,1721,1722,2174,2365,2489,2492,2682,2684 ^ 445 

635.762.953.1213.1598.2173.2300.2366.2876.3257.3642 ^ 443 

631.1199.1567.2159.2335.2796.2802.2972.2993.3738.3753 ^ 439 
607,1207,1591,2143,2351,2778,2780,2985,2988,3761,3762 ^ 431 

623.1183.1583.2167.2359.2794.2804.2969.2996.3737.3754 ^ 415 
829,1724,2490,2681,3131 ^ 254 
253,761,827,830,954,2297,2874,3133,3134,3260,3644 ^ 251 
247,799,823,2793,2801,2970,2994,3103,3119,3740,3756 ^ 239 
815,2777,2986,3127,3764 ^ 223 
575.2111.2745.2746.2748 ^ 191 



t = 8 

3783,3975 ^ 3015 

2967,3727 ^ 2791 

2775,2959,2983,3735,3751 ^ 2767 

2735,2743 ^ 2719 

2021,3534,3555,3917,3923 ^ 2006 

2019,3533,3557,3918,3926 ^ 2003 

1971,2541,2910,3501,3678 ^ 1779 

1949,2550,2933,3483,3691 ^ 1774 

1947,2549,2934,3485,3694 ^ 1771 

1943,2535,2903,3030,3045,3471,3663,3790,3811,3981,3987 ^ 1767 
1782,1965,1973,2526,2542,2909,2925,3499,3507,3675,3699 ^ 1758 
1963,2525,2926,3509,3702 ^ 1755 

1959,2511,2895,3021,3022,3495,3671,3795,3798,4003,4005 ^ 1751 
2493,2686 ^ 1723 

2479,2655,2782,2989,3763 ^ 1719 

1711,2463,2487,2671,2679,2798,2806,2973,2997,3739,3755 ^ 1695 
1532,1913,1916,3449,3450 ^ 1530 

1914,3452 1529 

1526.1885.1909.2012.2028.3419.3435.3546.3570.3945.3953 ^ 1518 
1523,1886,1907,2033,2034,3422,3437,3561,3564,3930,3932 ^ 1517 

1525.1883.1910.2010.2025.3421.3438.3548.3569.3948.3954 1515 
1901,2036,3443,3562,3929 1502 



continued on next page 



83 



t — 8 (continued) 
1501,1899,1902,2009,2026,3445,3446,3545,3572,3946,3956 ^ 1499 

1469,1659,1662,1786,1977,2429,2430,2556,2940,3513,3706 ^ 1467 
1463,1631,1655,2399,2415,3036,3052,3802,3826,4009,4017 ^ 1455 
1647,2423,3060,3818,3993 ^ 1439 

1470,1661,2427,3322,3897 ^ 1020 

1014,1773,1950,2547,2931,3307,3310,3486,3693,3867,3869 ^ 1013 
3309,3870 1011 
1005,1781,1974,2539,2907,3294,3315,3502,3677,3885,3891 ^ 990 
1006,1757,1966,2523,2923,3291,3318,3510,3701,3883,3893 ^ 989 
999,3027,3043,3279,3303,3789,3813,3855,3863,3982,3990 ^ 983 
957,1725,1726,2302,2491,2494,2683,2685,2877,3259,3643 ^ 766 

958.2299.2875.3262.3645 ^ 765 

955.2301.2878.3261.3646 763 

951.2287.2797.2803.2847.2974.2995.3247.3615.3742.3757 ^ 759 

927.2295.2795.2805.2871.2971.2998.3231.3631.3741.3758 ^ 751 
943,2271,2779,2781,2863,2987,2990,3255,3639,3765,3766 ^ 735 
2239,2623,2747,2749,2750 -> 703 
893,1278,1788,1853,1980,2554,2937,3195,3387,3514,3705 ^ 510 
894,1275,1785,1851,1978,2553,2938,3198,3389,3516,3708 ^ 509 
891,1017,1018,1277,1854,3197,3321,3324,3390,3898,3900 ^ 507 
887,1263,1823,3057,3058,3183,3359,3817,3820,3994,3996 ^ 503 
863,1271,1847,3034,3049,3167,3375,3804,3825,4012,4018 ^ 495 
879,1247,1839,3033,3050,3191,3383,3801,3828,4010,4020 -> 479 
639,1215,1599,2175,2367,2810,2812,3001,3004,3769,3770 ^ 447 
831,2809,3002,3135,3772 ^ 255 



t = 9 

3047,3791,3815,3983,3991 ^ 3031 

3799,4007 ^ 3023 

2807,2975,2999,3743,3759 ^ 2799 

2991,3767 ^ 2783 

3565,3934 ^ 2035 

2029,2037,2038,3550,3563,3566,3571,3931,3933,3949,3955 ^ 2014 

2030.3547.3574.3947.3957 ^ 2013 

2027.3549.3573.3950.3958 ^ 2011 
2023,3535,3559,3919,3927,4045,4046,4051,4054,4067,4069 ^ 2007 
1981,2558,2941,3515,3707 ^ 1790 
1979,2557,2942,3517,3710 ^ 1787 
1975,2543,2911,3038,3053,3503,3679,3806,3827,4013,4019 ^ 1783 
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t — 9 (continued) 
1951,2551,2935,3061,3062,3487,3695,3819,3822,3995,3997 ^ 1775 

1967.2527.2927.3037.3054.3511.3703.3803.3830.4011.4021 ^ 1759 
2495,2687,2814,3005,3771 ^ 1727 
1917,2044,3451,3578,3961 ^ 1534 
1533,1915,1918,2041,2042,3453,3454,3577,3580,3962,3964 ^ 1531 
1527,1887,1911,3423,3439,4058,4060,4073,4076,4081,4082 ^ 1519 
1903,3447,4057,4074,4084 ^ 1503 
1663,2431,3068,3834,4025 ^ 1471 
1022,1789,1982,2555,2939,3323,3326,3518,3709,3899,3901 ^ 1021 
3325,3902 ^ 1019 
3059,3311,3821,3871,3998 ^ 1015 

1007.3035.3051.3295.3319.3805.3829.3887.3895.4014.4022 ^ 991 
959,2303,2811,2813,2879,3003,3006,3263,3647,3773,3774 767 
895,1279,1855,3065,3066,3199,3391,3833,3836,4026,4028 ^ 511 



Appendix C 
Master Integrals 



The irreducible integrals that appear in the reduction of the three-loop quark form 
factor were given diagrammatically in fig. 3.1 (t=9), fig. 4.2 (t=7), fig. 4.1 (t=6) 
and fig. 5.2 (t=5). Here we shall collect results for these integrals as far as they are 
known in the literature. They are also hsted in [28, 29, 39]. 

For the integral Ag^i (ID=1790, see fig. 3.1 and eq. (3.7)), a result has recently 
been obtained [39]: 
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^46,1 (ID=1433) is integrable loop-by-loop. We may express the result using eqns. (5.5) 
and (5.6): 

A,l=H'-^r?)"'Vl>,2e,l. (C.4) 



^46,2 (ID=444) and ^6,3 (ID=429) are much more difficult to evaluate. The authors 
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of [29] found 
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The t — A and t = 5 topologies in fig. 5.2 read 
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(ID = 172) 
(ID = 158) 
(ID = 433) 
(ID = 412) 
(ID = 662) , 



with the /'s and g^s as defined in (5.5) and (5.6). 
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